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ABSTRACT 


NONLINEAR  OPTICAL  POLARIZABILITIES 
OF  MOLECULES 

by 

Charles  Kenneth  Miller 


Chairman:  John  F.  Ward 


The  technique  of  dc-electric-field-induced  second- 
harmonic  generation  (dcSHG)  has  been  used  to  investigate 
higher-order  polarizabilities  (hyperpolarizabilities)  of 
molecules  in  gases  and  vapors,  and  a model  has  been 
developed  to  describe  trends  in  dipole  moments,  polar- 
izabilities and  hyperpolarizabilities  for  series  of 
structurally  related  molecules. 

A review  is  given  of  the  relation  between  microscopic 
and  macroscopic  forms  of  electric  tensor  properties  of 
both  second  and  third  orders.  The  dcSHG  process  in  a 
focused  laser  beam  for  a particular  dc  field  geometry  is 
also  reviewed,  and  the  treatment  is  extended  to  include 


anomalously  dispersive  gases. 

The  experimental  apparatus  comprises  a Q-switched 
ruby  laser,  optical  components,  detectors,  and  signal 
processing  electronics,  as  well  as  a gas-handling  and  /* 

/f. 

pressure  measurement  system  which  permits  temperature-  /f 

A?  -S' 

dependent  gas-phase  observations  on  substances  which  ar 


D 

0 

liquids  at  room  temperature.  A review  is  made  of  the 
various  experiments  by  which  the  second-  and  third-order 
polarizabilities,  ratios  of  hyperpolarizability  components, 
and  wave-vector  mismatches  are  measured. 

Hyperpolarizabilities  and  hyperpolarizability  ratios 


t 


have  been  measured  for  the  molecules  H?,  N2,  0^,  CO,  NO, 
COg,  SFg,  NH^,  H2S,  H^O,  (CH^)20,  qnd  CH^OH.  The  results 
are  presented  and  compared  with  theoretical  and  experi- 
mental data  from  other  sources.  Agreement  with  most  other 
gas-phase  experimental  data  is  good,  and  where  this  is 
not  the  case,  we  believe  our  data  to  be  the  more  reliable. 
Special  points  of  interest  are  discussed  for  several  of 
the  molecules. 

The  ability  to  work  in  the  vapor  phase  with  sub- 
stances which  are  liquids  at  room  temperature  has  enabled 
the  acquisition  of  hyperpolarizability  data  on  several 
halogenated  methanes  which  are  liquids  at  room  temperature. 
Such  data  are  now  available  for  all  molecules  of  the  form 

CXnY4-n'  n = 1'Ut  X'Y  = F'C1'H* 

Correlations  of  electric  tensor  properties  among 

molecules  of  similar  structure  are  discussed.  The 

familiar  Bond  Additivity  Approximation  is  reviewed  and 

applied  to  the  dipole  moments,  linear  polarizabilities, 

and  hyperpolarizabilities  of  the  set  H^O^CHg^O,  and 


CH^OH,  and  to  the  set  of  halogenated  methanes.  A model 


of  molecular  electric  tensor  properties  is  developed  in 


ii 


si 


n 


D 

n 


.. . . 


r 


. . 

D 

0 


■ *- 


which  electrostatic  dipolar  inter-bond  interactions  are 
accounted  for  in  a self-consistent  manner.  A simplifica- 
tion of  the  model  is  applied  with  limited  success  to 
the  set  of  halogenated  methanes,  and  trends  in  the  sizes 
of  various  interaction  terms  are  noted. 
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CHAPTER  1 


INTRODUCTION 

The  response  of  an  isolated  molecule  to  a strong 
electric  field  is  of  considerable  theoretical  and 
experimental  interest.  The  ability  to  predict  linear 
and  nonlinear  polarizabilities,  which  are  measures  of 
this  response,  is  a stringent  test  of  molecular  orbital 
calculations,  and  accurate  data  are  necessary  for  compari- 
son. A knowledge  of  these  values  is  also  necessary  to  the 
experimentalist  who  may,  for  example,  wish  to  separate 
out  the  effects  of  interparticle  interactions  in  the 
polarizabilities  of  condensed  phases.  We  have  measured 
nonlinear  polarizabilities  o " a number  of  molecules  with 
the  aid  of  the  strong  optical  electric  field  of  a focused 
ruby  laser  beam,  and  we  have  developed  a model  which  is 
useful  in  correlating  the  permanent  dipole  moment  and 
linear  and  nonlinear  polarizabilities  in  groups  of 
similar  molecules. 

A molecular  system  subjected  to  a (time-dependent) 
electric  field  E possesses  an  electric  dipole  moment  p.(E). 
Disregarding  for  the  moment  the  tensorial  nature  of  the 
interaction,  the  relationship  can  be  written 


M-(E)  = u + XE 


(1.1) 


in  which  u is  the  permanent  dipole  moment  and  x is  the 
electric  dipole  polarizability,  x is,  in  general,  a 


1 


2 


function  of  the  applied  field,  and  if  it  is  analytic  for 
a range  of  E,  then  in  that  range  equation  1.1  may  be 
written 

U(E)  = u + aE  + X(2)E2  + X(3)E3  + ...  (1.2) 

where  a and  x^  are  independent  of  E (the  frequency 
dependence  of  a and  X^^  is  suppressed  here;  a more 
general  equation  will  be  presented  in  chapter  2).  The 
linear  polarizability  a is  responsible  for  the  familiar 
dielectric  properties  of  matter.  For  applied  fields  of 
optical  frequency,  a is  manifested  in  such  linear  optical 
effects  as  the  deviation  of  the  linear  index  of  refraction 
from  1 and  Rayleigh  scattering.  A complete  discussion 
of  these  phenomena  can  be  found  in  reference  1.  The 
quantities  x^»  n > 2,  are  the  nonlinear  polarizabilities, 
or  hyperpolarizabilities.  Non-zero  values  of  these 
coefficients  imply  that  the  induced  dipole  moment  has 
terms  proportional  to  second  and  higher  powers  of  the 
applied  field,  and  therefore  that  it  may  have  components 
at  frequencies  which  are  svuns  and  differences  of  the 
applied  field  frequencies.  Effects  due  to  these  higher- 
order  terms  in  the  induced  moment  fall  under  the  classi- 
fication of  non-linear  optics. 

Until  the  early  1960's,  available  optical  electric 
fields  of  any  particular  frequency  were  small  in  magnitude, 
and  effects  due  to  the  products  of  these  fields  were 
virtually  undetectable;  thus,  observed  ’’nonlinear  optical" 
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effects  were  limited  to  those  in  which  the  optical 
polarization  was  at  the  same  frequency  as  the  incident 
optical  field,  and  the  higher  order  terms  in  u(E)  contained 
products  of  the  optical  field  with  dc  fields.  Examples  of 

such  processes  are  the  (dc)  Kerr  effect  and  the  linear 

2 

electro-optic  effect.  With  the  advent  of  powerful 
lasers,  strong  electric  fields  become  available  at 
optical  frequencies,  permitting  the  observation  of  non- 
linear optical  phenomena  of  many  types. 3 

Of  interest  in  this  dissertation  is  dc-electric-f ield- 
induced  second-harmonic  generation  (dcSHG)  in  molecules  in 
the  gas  phase.  In  dcSHG,  the  x^E^  term  of  equation  1.2, 
in  which  the  field  has  components  at  optical  frequency  u> 
and  at  zero  frequency  (dc),  yields  a polarization  with  a 

component  at  the  frequency  2ao.  The  dcSHG  process  was  first 

4 

studied  in  gases  by  Mayer  _et  al.  Subsequently,  more 

5 6 

accurate  measurements  were  made  by  Finn  and  by  Bigio, 

and  the  experiments  described  herein  are  in  part  an 

extension  of  these.  We  have  expanded  studies  of  the  non- 

S 6 

linear  properties  of  a group  of  halogenated  methanes^* 
to  include  CCl^,  CHCly  CHjCl^,  CH^Cl,  CCl^F,  and  CCl^, 
so  that  data  are  now  available  for  all  twelve  molecules  of 
the  form  CX  Y^_  with  n = 1-4  and  X,Y  = F,C1,H.  In 
addition,  a number  of  other  small  molecules  were 
investigated:  H?,  Ng,  02,  CO,  NO,  C0?,  HgS,  H20,  (CH^O, 

CH^OH , and  SF^.  All  of  the  molecular  systems  studied 
here  are  simple  enough  to  be  potential  candidates  for 


4 


ab  initio  calculations  of  hyperpolarizabilities 
(calculations  already  exist  for  some),  and  data  from 
various  other  nonlinear  optical  experiments  are  available 
for  many  of  them.  In  addition,  several  of  the  molecules 
and  groups  of  molecules  have  features  of  particular 
interest,  to  be  discussed  in  due  course. 

In  our  experiments,  the  beam  from  a pulsed  ruby  laser 
is  focused  in  a gas  in  a region  of  dc  electric  field,  and 
the  second-harmonic  radiation  generated  at  the  focus  is 
detected  and  measured.  The  problems  of  relating  the 
nonlinear  susceptibility  of  a molecular  gas  to  the 
hyperpolarizability  tensors  of  the  molecules,  and  of 
determining  this  susceptibility  from  the  second  harmonic 
generated,  are  reviewed  in  chapter  2.  In  chapter  3 > 
attention  is  turned  to  the  experimental  apparatus,  and  a 
description  of  the  various  types  of  experiments  performed 
is  given. 

Experimental  data  and  derived  coefficients  for  the 
molecules  other  than  halogenated  methanes  are  presented  in 
chapter  4.  Comparisons  with  other  experimental  data  and 
with  available  theoretical  calculations  are  included,  as 
well  as  discussions  of  specific  points  of  interest  for 
several  of  the  molecules.  Chapter  5 contains  the  data, 
coefficients,  and  comparisons  with  experiments  and  theory 
for  the  halogenated  methanes. 

The  availability  of  the  data  for  all  of  the  halogenated 
methanes  facilitates  the  critical  evaluation  of  various 
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schemes  for  the  correlation  of  molecular  electric 
tensor  properties  among  homologous  molecules.  In  chapter 
6,  the  Bond  Additivity  Approximation,  a correlation  scheme 
in  which  inter-bond  interactions  are  ignored,  is  described 
with  examples  of  its  application  and  with  notice  taken  of 
its  drawbacks.  An  Interacting  Segment  Model,  in  which 
intramolecular  interactions  are  accounted  for  approximately, 
is  presented  and  applied  to  the  dipole  moments  and  to  the 
linear  and  second-order  polarizabilities  of  the  halogenated 
methanes.  The  results  are  compared  with  those  of  bond 
additivity. 

Chapter  7 contains  a summary  and  conclusions  concern- 
ing the  data  and  the  model. 


CHAPTER  2 
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THE  THEORY  OF  DCSHG  IN  MOLECULAR  GASES 
The  objective  of  the  dcSHG  experiments  described  in 
this  dissertation  is  to  gain  information  about  the  non- 
linear polarizabilities  of  isolated  molecules.  The 
experiments,  however,  must  actually  be  performed  on  an 
ensemble  of  molecules,  so  that  the  nonlinear  properties 
observed  will  be  ensemble  averages  of  those  individual 
molecules  and  will  therefore  in  general  have  different 
symmetry  properties.  In  particular,  a gas  is  macroscop- 
ically  isotropic,  even  though  an  individual  molecule  of  the 
gas  may  not  be.  In  addition,  dispersive  effects,  which  can 
be  ignored  for  a single  molecule  of  dimensions  much  shorter 
than  the  wavelength  o"  the  fundamental  or  second-harmonic 
radiation,  must  be  considered  when  the  light  traverses  a 
sample  of  gas  which  is  long  compared  with  the  wavelength. 

In  this  chapter,  the  dcSHG  process  is  reviewed  for  a 
microscopic  system,  after  which  the  relation  between  micro- 
scopic and  macroscopic  hyperpolarizabilities  in  a dipolar 
molecular  gas  is  examined.  The  chapter  is  concluded  with 
a brief  review  of  second-harmonic  generation  in  the  geometry 
of  these  experiments,  including  the  effects  of  dispersion 
and  a non-uniform  dc  electric  field. 

2. 1 Microscopic  Hyperpolarlzab  illties 

Equations  1.1  and  1.2  are  scalar  simplifications  of 
full  tensor  equations.  For  the  general  situation  of  applied 
field  and  induced  moment  with  several  frequency  components 
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it  is  useful  to  express  the  Fourier  amplitudes  of  the 
frequency  and  spatial  components  of  the  induced  moment  as 
functions  of  the  Fourier  amplitudes  of  the  frequency  and 
spatial  components  of  the  applied  field.  The  general 
equation  for  the  Fourier  amplitude  of  the  induced  dipole 
moment  per  molecule  at  frequency  oo  is 


•I. 


I K(-'-VaV 

. . +03  =05 

n rr 


f 


V 


05,  Ol 

X X ^ 4 t { ~ ^5_  J05,  , . • . >05  ) E . ...  K . 

J1l-**Tn  a 1 n ll  1n 


(2.1) 


where  summation  over  repeated  indices  is  implied.  The 

05|j, 

Ej  are  Fourier  amplitudes  of  the  time-dependent 
components  of  the  applied  electric  field  with  frequency 
o>m  and  polarization  vector  j,  given  by 

e/V)  - Ke/””  e'*®”1  + Ej'®"1  e1®”1)  (2.2) 

where  E“a'  is  the  complex  conjugate  of  Ef°.  The  definition 
of  K(-o5fT;o51> . . . ,05n)  is  discussed  by  Orr  and  Ward. 
Basically,  K includes  factors  from  Fourier  decomposition 
of  fields  and  from  distinguishable  permutations  of  field 
frequency  labels,  with  the  results  that: 

1.  A single  ordering  of  field  frequencies 
may  be  used,  i.e.,  one. need  only  write. 


...  . 
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say,  x(-0Jcr;a)1,a^2)  and  not  x(-o> 
in  addition; 

2*  lim  X( • • • »^n)=X  (0;0, . . . 0) ; and 

, . . . ,£Dn-*0 

3.  K(0;0, . . . ,0)=1  for  all  orders  of  hyper- 

polarizability. 

For  a molecule  subjected  to  an  optical  electric  field 
JSa  and  a dc  field  E°,  the  dipole  moment  amplitude  at 
frequency  2o>  to  third  order  in  the  applied  fields  is, 
from  equation  2.1, 


2od  1 


0) 


H = 2 Xikje(-2cD;cD,CD)Ek  E£ 


(2.3) 


where  the  factors  1/2  and  3/2  are  the  values  of 
K( -2a>;a>,cu)  and  K( -2a.'; 0,co,cjd)  respectively.  To  get  a rough 
idea  of  the  relative  sizes  of  the  contributions  to 
it  is  useful  to  examine  the  form  of  time-dependent  quantum- 

7 

theoretical  perturbation  expansions  for  the  non-resonant 
second-  and  third-order  contributions  to  equation  2.3. 

These  expansions  contain  terms  like  the  following: 


(second  order  term)  = 


<XlWVmn<Vna  c ^ , 


h (wm  +2a)  (uu  -u ) 
v mg  > ' ng  ' 


(2.4) 


I 
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u 


Li 


(third  order  term)  = 


-»e4  <Xi>gD<X1>P.,<XR>mn<Xt>ng 

(V+2',')(u”i8‘2,')(Un6'l") 


E 


J 


u 


(2.5) 


U 


i 


n 


where  the  co  are  molecular  transition  frequencies  and 
rs 

e<X  > are  dipole  matrix  elements  between  states, 
q rs 

To  obtain  an  order-of-magnitude  estimate  of  the  relative 
sizes  of  these  terms,  we  first  replace  the  matrix  elements 
by  a typical  molecular  dimension  r.  Further,  we  assume 
that  the  as  are  all  much  larger  than  to  or  2o>,  and  can  all 

i o 

be  replaced  by  one  effective  frequency  a'Q.  Then,  noting 

2 

that  the  electron  energy  ha>o  is  approximately  e /r  and 

p 

the  typical  intramolecular  field  Emol  is  of  the  order  e/r“, 
we  see  that  the  third-order  polarization  terms  differ  from 
the  second-order  by  a factor  of  E°/Emo^.  With  a typical 

o 7 

r of  1 A,  Emol  is  of  the  order  of  10  esu;  thus  u 
(second-order)  is  much  larger  than  q2a'( third-order ) for 
electric  fields  E°  which  can  be  sustained  in  gases  at 

p 

ordinary  pressures  without  breakdown  (~  10"  esu).  As 
will  be  seen  in  the  next  section,  however,  incomplete 
molecular  alignment  causes  the  polarizations  arising  from 
second-  and  third-order  molecular  SHG  processes  that  are 
actually  observed  in  the  gas  phase  to  be  of  roughly  the 
same  magnitude. 


10 


2.2  Relating  Microscopic  and  Macroscopic  Hyperpolarlza- 

bilitles 

In  an  unperturbed  molecular  gas,  the  orientation  of 
the  coordinate  frame  of  a given  molecule  with  respect  to 
the  laboratory  coordinate  frame  is  completely  random.  In 
what  follows,  the  alignment  of  molecules  by  an  electric 
field  is  treated  classically.  This  approach  is  justified 
because  the  thermal  energy  kT  at  room  temperature  is 
sufficiently  larger  than  the  rotational  constants  for  all 
dipolar  molecules  studied  that  the  classical  limit  of 
orientational  averaging  is  consistent  with  a quantum 
mechanical  rotational  state  distribution.  If 
the  gas  is  subjected  to  a uniform  electric  field  E and 
transient  effects  are  allowed  to  die  out,  distribution  of 
molecular  orientations  will  contain  the  Boltzmann  weighting 
factor 

g = exp(-uAT),  (2.6) 

where 

1 

u = - < r u(?E).EdO  , (2.7) 

o 

k is  Boltzmann’s  constant,  and  the  angular  brackets 
indicate  an  average  over  a time  long  compared  with  the 
orientational  relaxation  time  of  the  molecule.  The  dipole 
moment  will,  in  general,  contain  contributions  from  all 
orders  of  polarizability.  The  magnitude  of  u is,  to  a 
good  approximation.  Just  u.E°,  where  u is  the  magnitude 
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*|  o 

of  the  permanent  dipole  moment  of  the  molecule  (~  10~  esu) 

and  E°  is  the  magnitude  of  the  dc  component  of  E.  With  an 

E°  of  ~ 40  esu,  the  orientational  energy  is  about  4xl0~^erg, 

i -i4 

compared  to  kT  of  4x10  erg  at  room  temperature.  Thus 
we  may  write 

g ~ 1 - u/kT.  (2.8) 

Q 

The  orientational  average  of  j_l  03  is  given  by 


, 2cov  ^ 

<H  > = x 


g do 


(2.9) 


where  integration  over  Q indicates  formally  the  integration 
over  all  orientations  of  molecular  coordinates  with  respect 
to  laboratory  coordinates.  In  equation  2.9  M-  has  the 
form 

2o)  1 , / 0 x . x . r,  OOt-,  0) 

^i  “ 2 Xik£  ("2a';a5,a^$Fi*Hk$M£EH  EM 


+ I Xijk£  (“2a';0'“'^*Fi$Gj$Hk$M£EG  EH  EM  (2-10) 


where  i,j,k,£  refer  to  the  molecular  coordinate  frame, 
F,G,H,M  refer  to  the  laboratory  frame,  and  the  $Nr  are 
direction  cosines  between  the  molecular  and  lab  frame 
coordinate  axes.  The  only  terms  retained  in  equation  2.9 
are  those  of  up  to  first  order  in  E°  and  up  to  second 
order  in  e5°.  The  only  term  in  u contributing  to  these 
coefficients  is  m._.$g.jEg0  (another  term  of  comparable  size 
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in  this  experiment  would  be  (1/2)  xki  (-cdjcdH^^E^E^  ,but 
it  does  not  survive  the  orientational  averaging).  The 
integral  in  equation  2.9  has  been  performed^  and  yields 
the  result 


..  2u>  _ 3/  (3)  . 1 . v (2 ) \E  o cd„  cjo 

UF  ~ 5\XFGHM  9kT  4XFGHM/EG  EH  ^ ' 


(2.11) 


where 


X FGHM-  ^Fi$G,1$Hk#M£>xiJk£^‘2a3;0,a'’,f°) 


(2.12) 


^FGHM"  <$Fi$Gz$Hk$Mi>M-zxikt^'2cu;a3,a^  (2-13) 

The  are  averaSes  of  products  of  direction 

cosines  over  all  molecular  orientations,  which  have  been 

8 A 

tabulated  by  Cyvin  et  al. , and  z is  the  direction  of  the 
molecular  dipole  moment,  as  discussed  in  section  3.2c. 

It  should  be  noted  that  all  contributions  to  u2cD 
of  second  order  in  the  electric  field  average  to  zero  in 


equation  2.9>  so  that  one  may  write 


..  2cd  _ 3 e p Op  uv,  cd 
UF  “ ? xFGHMEG  EH  EM 


(2.14) 


where  the  effective  hyperpolarizability  per  molecule,  x > 


is  given  by 


e (3)  u (2) 

XFGHM  = XFGHM  §kT  XFGHM 


(2.15) 
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The  isotropy  of  the  gas  phase  places  strong  restrictions 
on  Xe.  For  example,  given  the  coordinate  transformation 
X'=-X,  isotropy  demands  that 

ve  - ve  (2  16) 

XX'YYY  ~ XXYYY  ' * ' 

Now  considering  the  process 


2cd, 

X 


_ 3 6 pi  Op,  COp,  CD 

- ? xxyyyey  ey  ey 


(2.17) 


the  electric  fields  in  the  ^ direction  do  not  change  sign 

under  the  coordinate  transformation,  but  the  dipole 

moment,  having  the  parity  of  an  electric  field  in  the 
A 

X direction,  does  change  sign,  so  that 


2cd 


T ^ “ 9 X Y» 

3 e 


2a\  3 „ e E o cdf  cd 

X'YYY  Y rjY  Y 


- S *XYYYEY°W  - 


2a)\ 


(2.18) 


It  follows  that  <|ix  >,  and  therefore  XvyyY’  are  7ero • 

By  considering  similar  coordinate  transformations  and 
exchange  symmetry  of  the  two  electric  fields  at  frequency 
co,  it  follows"^  that  there  are  only  two  independent  non- 
zero elements  of  XpGHNT  and  therefore  of  and  ^FGHM' 

and  that  they  can  be  chosen  to  be  FGHM=YYYY  (to  be  denoted 
as  || ) and  FGHM=YYXX  (to  be  denoted  i).  The  symbols  n and 
i thus  stand  for  the  situations  in  which  the  applied 
optical  electric  field  polarization  is  respectively 
parallel  to  and  perpendicular  to  the  dc  field  direction. 

X^  and  x(2)  can  be  evaluated^  in  terms  of  molecular-frame 
' ii  > A ii  > 
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components  using  equations  2.12  and  2.13  to  give 

x!,3)=  T5  + V5  ^aar)^-2^0 t2'19) 

xl3)=  T5  x<jaTn(-a“;0'a>'aj)  “ T5  XaW -2to;°^,co)  (2.20) 

xIi2)=  I XTnzC-2^^)  + | (2-21) 

xf)=  | ’W"2®*®'®)  (2*22) 


(2.22) 


where  summation  over  repeated  indices  is  implied. 

The  nonlinear  polarization  of  the  medium  is  the 
average  induced  nonlinear  dipole  moment  of  each  molecule 
times  the  number  density  of  the  molecules: 


p 2o)  _ r e F Op  ODp  05 

PF  - pl*FGhmeg  eh  ^ 


(2.23) 


where  L is  Loschmidt's  number  (2.69xl019  cm-^)  and  p is 
the  number  density  of  the  gas  in  units  of  Loschmidt's 
number  (p  is  a dimensionless  quantity  and  is  equal  to 
1 for  an  ideal  gas  at  STP).  In  equation  2.23,  no 
distinction  is  made  between  the  applied  electric  field 
and  the  local  field.  The  inequivalence  of  these  fields 
is  important  in  condensed  phases;9  however,  in  gases  at 
the  pressures  used  in  these  experiments,  intermolecular 

p. 

interactions  are  small,  and  the  local  field  can  be 
assumed  to  be  equal  to  the  applied  field. 

2 . 3 dcSHG  in  a Focused  Gaussian  Beam 

The  polarization  at  the  optical  frequency  2o>  given 


»■  4 **  * .» 


• « • 
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by  equation  2.23  will  emit  electromagnetic  radiation  at 

that  frequency.  The  pattern  and  intensity  of  this 

emission  depends  on  the  spatial  distribution  of  E° 

and  pf*5,  and  of  the  dispersion  of  the  gas.  Dependence  of 

dcSHG  on  the  parameters  of  this  experiment  is  discussed  by 
c,  10 

Finn  and  by  Bigio  et  al.  , and  only  the  major  points  are 
sketched  here. 

The  propagation  of  the  second-harmonic  beam  is  des- 
cribed by  the  inhomogeneous  wave  equation 

P 2 

_ _ 1-205/  x . W t-2<jd/v  i,  4q>  Dn£  , v 
v x v x F (r)  + e — E (r)  = -%  — P (r) 

c c 

(2.24) 

where  P (_r1  is  the  driven  second-harmonic  polarization 
given  by  equation  2.23.  The  incident  optical  field  for  a 

A 

lowest-order  Gaussian  beam  propagating  in  the  z direction 
and  polarized  in  the  y direction  is  given  by 


„ a Qika'b§/2  / >o) / 2,  2. . 

Ef“(r)  . vFf  5 — exp  (-K  +y-l)  ' 

1+15  b(l+i5) 


(2.25) 


where  b is  the  confocal  parameter,  k^  is  the  wave  vector 
at  frequency  oo,  and  the  dimensionless  parameter  5 , equal 
to  2z/b,  is  introduced  for  convenience.  The  second- 
harmonic  field  is  assumed  also  to  be  a Gaussian  beam 
coaxial  and  confocal  with  the  incident  beam,  with  the 
same  b and  with  k^  replaced  by  k' ,J),  but  having  an  addi- 
tional factor  of  a slowly-varying  z dependence,  reflecting 
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the  (assumed)  slow,  uniform  generation  of  second 

A 

harmonic  as  the  beam  proceeds  in  the  z direction.  All 
that  is  needed  to  complete  the  definition  of  equation 
2.25  is  the  specification  of  the  dc  field. 

In  these  experiments,  the  dc  electric  field  is 
generated  by  maintaining  a voltage  difference  between 
two  long,  parallel  cylindrical  electrodes  each  of  radius 
a,  with  a distance  / between  centers,  as  shown  in  figure  1. 
The  electric  field  on  the  z axis  is  calculated^  to  be 


E°(z)  = $ 


d cosh"1 (//2a)  (1  + 4z  /d  ) 


(2.26) 


where 


(/2  - wy/2 


(2.27^ 


is  the  distance  between  equivalent  line  charges  and  Vq 

is  the  potential  difference  between  the  electrodes.  The 

positional  relationship  of  the  dc  electric  field  with 

respect  to  the  incident  beam  is  shown  in  figure  1.  The 

beam  is  perpendicular  to  the  plane  of  the  electrode 

centers  and  equidistant  from  them,  with  the  waist 

of  the  beam  located  a distance  z„  from  the  electrode 

o 

plane.  Since  the  beam  is  of  small  extent  in  the  x and  y 
directions  near  the  electrodes,  it  is  a good  approximation 
to  say  that  E°(r)  = E°(z). 

Using  the  aforementioned  approximations,  and  adding 
those  of  small  far-field  diffraction  angle,  negligible 
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fundamental-beam  depletion,  and  path-length  in  the  gas 
much  longer  than  b on  either  side  of  the  focus, 
equations  2.24-2.26  yield11  the  ratio  of 
generated  second-harmonic  power  to  the  square  of  the 


incident-beam  power: 


^7-  TgAA?  |^1|2  (E°(0))g  r(p)B(t)Z(zo)  (2.28) 


Here,  AkQ  is  the  wave-vector  mismatch  defined  by 

4k  . 2k  “ - k 2<“ 
o o o 


(2.29) 


where  the  subscript  o indicates  that  a quantity  is 
evaluated  at  o = 1.  The  dimensionless  functions  in 


equation  2.28  are  as  follows: 


r(p)  . [f  exp  (^i)]2  , 


B(b)  = 


Z(zc)  = 


- + / - 

d v b 


(2.30a) 


(2.30b) 


(2.30c) 


1 + (-7T2) 


The  quantity  pQ  in  equation  2.30a  is  equal  to  2/d|Ako|. 
Maximum  values  of  1 are  attained  for  the  functions  in 
equations  2.30a,  b,  and  c for  p equal  to  pQ,  b equal 
to  d,  and  zq  equal  to  zero,  respectively.  When  these 





However,  the  experimental  signal  depends  on  such  things 
as  alignment,  photomultiplier  quantum  efficiency,  and  so 
on,  which  are  susceptible  to  slow  random  drifts.  These 
factors  will  be  discussed  further  in  section  3.3. 

We  have  seen  that  a molecular  gas  subjected  to  a 

dc  electric  field  and  irradiated  by  an  optical  field  can 

have  a polarization  at  twice  the  optical  frequency  as  a 

consequence  of  the  second-  and  third-order  nonlinear 

(?)  ( 3^ 

polarizabilities  xv  and  xv  ' of  the  molecules.  Further, 
we  have  seen  how  this  polarization  will  radiate  at  the 
second  harmonic  frequency  for  a particular  configuration 
of  dc  and  optical  fields.  This  effect  can  be  exploited 
in  measuring  the  components  of  the  molecular  hyperpolariza- 
bilities. The  details  of  such  experiments  are  examined  in 
the  next  chapter. 


CHAPTER  3 

THE  DCSHG  EXPERIMENTS 

The  relation  embodied  in  equation  2.31  indicates 
how  measurements  of  dcSHG  in  molecular  gases  can  be  used 
to  determine  the  nonlinear  polarizabilities  of  the  molecule. 
In  the  first  section  of  this  chapter,  the  experimental 
apparatus  for  making  such  measurements  is  described.  In 
the  second  section,  the  experiments  performed  with  this 
apparatus  are  explained  individually.  The  third  section 
contains  a discussion  of  the  data  and  data  analysis,  in- 
cluding an  indication  of  statistical  error  limits  and 
comments  on  possible  systematic  errors. 

3. 1 The  Experimental  Apparatus 

The  apparatus  for  the  dcSHG  experiments  is  similar 
to  that  described  in  detail  by  Bigio,^  and  is  illustrated 
in  figure  2.  A brief  sketch  is  given  here  of  the  various 
components,  with  differences  from  Bigio’s  apparatus  noted. 

The  fundamental  optical  electric  field  is  generated 
by  a flashlamp-pumped  actively  ©-switched  ruby  laser,  Korad 
model  K-15,  operating  at  69^3  The  laser  cavity  contains 
a 1.5  mm  aperture  to  reduce  higher-order  transverse  modes, 
simultaneously  limiting  the  power  to  1-2  MW  in  a pulse  of 
~ 20  ns  duration.  The  laser  is  fitted  at  its  output  with  a 
Glan-Laser  polarizer,  manufactured  by  the  Karl  Lambrecht 
Corporation,  to  ensure  accurate  linear  polarization. 

20 
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Figure  2.  Schematic  Diagram  of  Experimental  Apparatus. 
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Approximately  10$  of  the  fundamental  beam  energy  is 

reflected  by  a pellicle  into  the  monitor  box.  Here,  a 

small  portion  of  this  reflected  monitor  beam  is  itself 

reflected  from  the  surface  of  a red  filter  onto  the  face 

of  a photodiode,  providing  a trigger  pulse  for  data- 

acquisition  electronics.  The  remainder  of  the  monitor 

beam  passes  through  a quartz  crystal  plate,  where  it 

( 

generates  second  harmonic  via  the  xv  ' of  the  crystal.  The 
fundamental  frequency  is  filtered  out  of  the  monitor  beam 
by  passage  through  5 cm  of  a saturated  solution  of  copper 
sulfate  in  water,  and  the  second  harmonic  is  further 
isolated  by  a narrow-band  interference  filter  with  peak 
transmission  centered  at  3^71  A.  The  second  harmonic  beam 
is  detected  by  a photomultiplier  tube  with  a dynode  circuit 
optimized  for  linear  operation  with  short,  high-current, 
pulses.  The  output  of  the  photomultiplier,  which  we  will 
call  the  monitor  signal,  is  proportional  to  the  square  of 
the  instantaneous  power  of  the  laser  beam. 

The  remaining  90$  of  the  beam  passes  through  a half- 
wave plate  if  rotation  of  the  plane  of  polarization  is 
required.  The  beam  then  passes  through  a 26  cm  focal 
length  lens,  followed  by  a red  filter  to  eliminate  UV  light 
from  the  laser  flashlamp  and  from  possible  incandescence 
at  the  surfaces  of  optical  components,  and  is  focused  in 
the  experimental  cell. 

The  experimental  cell  is  stainless  steel,  with 
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anti-reflection-coated  fused  silica  windows  sealed  to  it  by 
means  of  o-rings.  Criteria  for  o-ring  material  selection 
are  resistance  to  high  temperature  and  compatibility  with 
the  particular  gas  under  investigation.  Among  materials 
used  were  silicone  rubber,  ethylene  propylene.  Teflon,  and 
Viton.  The  cell  houses  the  electrodes  described  in 
section  2.3j  consisting  of  2-inch  long  polished 
stainless  steel  rods,  1/16  inch  in  diameter  and  approximate- 
ly 1/4  inch  apart.  The  cell  is  wrapped  with  heating  tapes, 
and  its  temperature  is  controlled  to  within  1°C  by  a 
proportioning  controller.  It  can  be  filled  with  gas  or 
vapor  or  evacuated  through  the  gas-handling  system  described 
below. 

The  laser  beam,  along  with  the  second  harmonic 
generated  in  the  interaction  region  of  the  cell,  proceeds 
to  the  detector  box,  in  which  the  fundamental  is  filtered 
out  and  the  second  harmonic  is  detected  in  the  same  manner 
as  in  the  monitor  box.  The  output  of  the  detector  box 
photomultiplier  will  be  referred  to  as  the  harmonic  signal. 
The  harmonic  and  monitor  signals  are  amplified  as  required 
and  are  gated  by  the  trigger  pulse  to  separate  integrators. 
The  integrated  charges  are  digitized,  and  the  digital 
values  are  punched  on  paper  tape  for  later  analysis  by 
computer. 

The  gas  handling  system  has  been  modified  from  that  of 
Bigio  to  allow  gas-phase  measurements  on  substances  which 
are  liquids  at  room  temperature.  The  room-temperature  and 
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high-temperature  gas  handling  systems  are  illustrated 
schematically  in  figure  3-  The  room-temperature  system 
includes  a Bourdon  tube  pressure  gauge  and  mercury  mano- 
meter for  pressure  measurement,  and  uses  brass  valves  with 
plastic  packing,  with  copper  tubing  for  interconnection. 

In  the  high-temperature  section,  a high-temperature 
capacitance  manometer  is  used  to  measure  pressures, 
interconnecting  tubing  is  made  of  stainless  steel  and  is 
heated  to  temperatures  above  those  used  in  the  cell  to 
ensure  against  condensation,  and  all-stainless-steel  bellows 
valves  with  heated  bonnets  are  used.  The  liquid  sample  is 
kept  in  a separately  heated  chamber  at  a temperature  such 
that  the  vapor  pressure  is  above  the  highest  pressure  to 
be  used  in  the  experiment.  Waste  vapors  are  condensed  in  a 
separate  cooled  chamber. 

The  capacitance  manometer  was  constructed  for  these 
experiments  to  measure  pressures  in  the  range  0-2  atmos- 
pheres absolute,  at  temperatures  in  excess  of  300°C.  A 
cutaway  drawing  of  the  high-temperature  pressure  sensing 
element  and  a schematic  diagram  of  the  measuring  circuit 
appear  in  figure  4.  The  pressure-dependent  variable  is 
the  capacitance  between  a stainless  steel  diaphragm  and  a 
fixed  metal  plate.  The  diaphragm  forms  one  wall  of  a 
capsule,  the  interior  of  which  is  connected  to  the  vapor 
being  measured.  The  other  side  of  the  diaphragm  is  at 
atmospheric  pressure.  The  thickness  of  the  insulator 
between  the  capsule  and  the  fixed  plate  was  chosen  so 
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Figure  3.  Room-Temperature  and  High-Temperature 
Clas-Handling  Systems. 


27 


that  the  coefficient  of  thermal  expansion  of  the  gap  would 
be  zero.  While  this  goal  was  not  achieved  exactly,  control 
of  the  sensor  temperature  to  1°C  was  sufficient  to  limit 
thermal  drift  to  an  acceptable  level.  The  capacitance  of 
the  sensor  head  forms  one  leg  of  an  ac  bridge  circuit;  the 
rest  of  the  bridge  consists  of  a fixed  silver  mica  capacitor 
and  the  two  legs  of  a Singer  Ratiotran  ratio  transformer, 
which  is  effectively  an  inductor  with  a finely  variable 
tap.  The  bridge  is  driven  with  a 10  kHz  sine  wave  from  a 
low-distortion  Wien  bridge  oscillator.  Any  out-of-null 
signal  from  the  bridge  is  detected  with  a simple  lock-in 
amplifier  circuit.  If  y is  the  ratio  of  the  inductance  of 
the  upper  leg  of  the  transformer  to  that  of  the  lower  leg, 
then  when  the  bridge  is  nulled,  y is  proportional  to  C , 
the  capacitance  of  the  sensor.  This  capacitance  unfor- 
tunately does  not  vary  linearly  with  the  pressure  difference 
across  the  diaphragm.  However,  one  can  expand  the  pressure 
difference  p as  a power  series  in  y 

P = a1(y-yQ)  + a2(y-yQ)2  + a3(y-yQ)3  , (3.1) 

and  likevise  y as  a function  of  p 

(y-yQ)  = b1p  + b2p?  + b3p3  (3.2) 

where  yQ  is  the  value  of  y at  p=o,  i.e.,  when  the  vapor 
is  at  atmospheric  pressure.  (Equations  3-1  and  3.2  are 
used  for  respectively  measuring  and  setting  pressure.) 
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Coefficients  and  can  be  chosen  such  that  these 
functions  fit  the  actual  p vs.  y data  to  better  than  2 mm  Hg. 
Drift  and  hysteresis  amount  to  no  more  than  5 mm  Hg.  for 
absolute  pressures  up  to  2 atm,  permitting  measurements  of 
|AkQ|,  as  described  below,  to  better  than  1 %.  The  sensor 
element  is  kept  continuously  at  a temperature  higher  than 
the  maximum  temperature  to  be  used  in  the  experimental  cell. 
If  the  sensor  temperature  is  changed,  a new  calibration  of 
the  a^  and  b^  is  made  before  the  gauge  is  used  again. 

3. 2 Experiment  Descriptions 

Ignoring  statistical  and  systematic  fluctuations 
(discussed  in  section  3 of  this  chapter)  the  harmonic 
signal  is  proportional  to  I2a>  and  the  monitor  signal  is 
proportional  to  (F1)2.  We  call  the  ratio  of  the  harmonic 
to  the  monitor  signals  simply  the  signal,  denoted  by  S. 

Given  the  interaction  geometry  of  figure  1,  we  see  from 
equation  2.31  that 


S = C 


IsfeJ  V 


(33) 


with  the  density  set  to  PQ  as  discussed  in  section  2.3. 

The  following  experiments  exploit  this  relation  to  determine 

X(3]  and  X{2]. 

\\>l  II » A 

3.2a  lx*l 

In  this  experiment  measurements  of  S for  a particular 


gas  are  made  alternately  with  measurements  of  S for  a 


reference  gas,  all  at  a particular  temperature.  Measure- 
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ments  relative  to  a reference  gas  are  favored  over  absolute 

measurements  because  of  difficulties  in  calibrating  C in 

equation  3-3  absolutely,  as  discussed  in  section  3 of  this 

chapter.  In  a typical  data  run  at  a fixed  temperature,  5 

subruns  of  40  laser  shots  each,  with  the  experimental  gas 

in  the  cell,  are  interleaved  with  5 sub runs  in  which  the 

reference  gas  is  in  the  cell.  Helium  was  used  as  the 

12 

reference  gas  bv  Ward  and  Bigio,  but  in  their  experiment, 
C was  calibrated  only  at  room  temperature  in  a cell  which 
could  withstand  the  high  pressure  at  optimum  density  for 
helium  (approximately  600  psi).  However,  for  substances 
which  are  liauids  at  room  temperatures,  C had  to  be 
calibrated  at  elevated  temperatures,  and  for  some  gases 
investigated,  the  change  in  |xe|  over  the  accessible 
temperature  range  was  comparable  to  systematic  drifts 
over  the  time  required  to  make  the  measurements,  so  it  v as 
found  desirable  to  calibrate  C at  various  temperatures.  For 
these  situations  S had  to  be  measured  for  the  reference  gas 
in  the  more  fragile  heated  cell,  which  contains  delicate 
ceramic-to-metal  seals,  so  a gas  with  a lower  pressure  at 
optimum  density  was  used.  Methane  was  chosen  as  the 
alternate  reference  because  of  its  ready  availability  and 
reasonable  (~  1 atm)  pressure  at  optimum  density.  Methane 
was  compared  carefully  with  helium  (see  section  4.2)  in  a 
high-pressure  cell.  It  should  be  noted  that  neither  helium 
nor  methane  has  a permanent  dipole  moment,  so  that  from 
equation  2.17  each  has  a x6  which  is  independent  of 
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temperature. 

Using  equation  3.3*  with  primes  to  indicate  reference 
gas  values,  we  can  write 


\1/2 

Ak 

o 

) 

"EE  r 
o 

(3.^) 


An  experiment  described  below  is  used  to  determine  AkQ  for 
the  experimental  gases  and  for  methane.  For  helium,  AkQ 
can  be  determined  from  accurate  data  on  the  dispersion  of 

1"5 

the  index  of  refraction: 


Ako(He)  = 0.09^6  cm' 


(3.5) 


To  establish  an  absolute  scale  for  the  magnitude  of 

0 1 ii 

X , we  use  a value  calculated  by  Sitz  and  Yaris  which  is 

thought  to  be  accurate  to  1$: 


X6  = X„„~„(“2,-ui0,a>,a))  = +3. 79x10" ^esu/helium  atom 

||  ZZZa 

(3.6) 

3.2b  Sign  of  xe 

According  to  equation  3.3>  the  signal  is  independent 

6 6 
of  the  sign  of  x • However,  the  relative  signs  of  x for 

two  different  molecules  can  be  established  by  measuring 

signals  generated  in  binary  mixtures  of  gases.  J If  the 

X6  are  of  opposite  sign,  cancellation  will  lead  to  zero 

signal  for  some  ratio  of  gas  densities.  In  practice,  the 

signal  is  measured  with  a fixed  density  of  reference  gas, 

and  the  density  of  experimental  gas  is  varied  through  a 


p . 

.*'*•.*  * , '' 
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range  including  the  expected  cancellation  value.  The 

variation  of  signal  with  density  of  experimental  gas  can 
5 

be  predicted  for  both  cancellation  and  no  cancellation, 
and  is  quite  different  for  the  two  cases;  in  addition,  as 
is  shown  in  the  Appendix,  the  cases  with  AkQ>0  (of  which 
none  were  encountered  in  these  experiments)  would  exhibit 
variations  quite  different  from  those  of  normally  dispersive 
gases. 

Two  possible  references  for  absolute  determination  of 

the  sign  of  xG  may  be  employed.  Either  the  (positive)  sign 

14 

of  the  calculated  value  for  helium  may  be  used,  or  one 

may  assume  that  the  experimentally  determined  (positive) 

16 

sign  of  the  Kerr  coefficient  for  argon-  is  the  same  as 
that  of  the  dcSHG  coefficient.  Binary  mixture  experiments1^ 
have  verified  that  the  dcSHG  coefficients  of  helium  and 
argon  are  of  the  same  sign,  so  that  there  is  no  inconsist- 
ency in  these  methods. 

A least-squares  fit  of  the  absolute,  signed  values  of 
X as  a linear  function  of  1/T  may  now  be  made  to  determine 
the  parameters  in  the  u version  of  equation  2.15,  which  can 
be  written  as 


+ 


M- 

WT 


The  best-fit  line  yields  the  values 


X 


(3) 

it 


lim  xe 


1/T-+0 


(3.7) 


* 


(3.8) 
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(2)  _ 9k  dx' 


(3.9) 


We  see  from  equation  3-7  that  since  the  sign  of  x (T) 

( o) 

has  been  measured,  the  sign  of  |XXV  ' is  determined 

(2) 

unambiguously.  We  choose  to  assign  to  x the  sign 
(2)  ’ 

measured  for  4XV  . ThiS'-choice  is  the  equivalent  of 
choosing  the  molecular  z-axis  to  have  the  same  direction 

and  sense  as  the  permanent  dipole  moment;  that  is,  the 

— + 

molecular  z-direction  is  defined  as  the  ► direction  of 

the  dipole  moment.  Hence,  4 is  equal  to  4 , which  is  always 

z 

positive.  For  some  of  the  molecules,  further  detailed 
discussion  of  results  requires  a knowledge  of  the  direction 
and  sense  of  the  dipole  moment  in  relation  to  molecular 
geometry.  These  cases  will  be  examined  in  chapter  6. 

3.2c  xx 

/"s')  (o} 

As  noted  in  section  2.2,  xK  ’ and  xv  each  have  two 

independent  nonzero  components,  denoted  by  u and  x.  Rather 

( -a\  (2) 

than  express  our  results  in  terms  of  xx  and  xx  directly, 
we  have  found  it  convenient  to  deal  with  the  n components 
and  the  ratios  R^  and  R^\  which  are  defined  as 
X<2>*<2>  and  x|^  Vx{^  > respectively.  If  there  is 
complete  exchange  symmetry  among  the  subscripts  of  the 


molecular  x ^ j ^ or  then  . > n , ^ , 

will  be  equal  to  3*  as  can  be  seen  by  examining  equations 

17 

2.19  - 2.22.  This  has  been  shown  to  be  the  case  to  within 
experimental  uncertainty  for  the  inert  gases.  Thus,  it  is 
the  deviations  of  R^  and  R^  from  3 which  will  be  of 


R^2)  or  R^3) 


or  RVJ',  respectively. 
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interest  experimentally. 

(p)  (-3) 

Two  methods  were  used  to  determine  Rv  and  Rw;. 

0 

Both  methods  involve  measuring  the  quantity  R (T)  defined 

by 


Re(T)  = X®(T)/X®(T)  . 


(3.10) 


(experimental  determination  of  this  quantity  involves 

relative  measurements  of  the  signal  as  a function  of  the 

angle  between  the  beam  polarization  direction,  as  adjusted 

by  the  half-wave  plate  in  figure  2,  and  the  direction  of 

the  dc  field,  as  described  in  detail  by  Bigio  and 

Ward  ).  The  first  method  is  direct,  and  depends  on  the 
0 

measurement  of  xx  at  various  temperatures.  The  analogs  of 
equations  3.8  and  3.9  are  then  used  to  determine  and 

after  which  R^^  and  R^^  are  calculated  directly. 

In  order  to  get  an  accurate  absolute  scale  for  x®>  the 
ratio  R (T)  is  measured  at  a single  temperature.  This 
procedure,  although  time-consuming,  is  straightforward,  and 
was  actually  used  for  the  first  gases  measured. 

In  the  second  method,  which  was  used  on  the  halogenated 
methanes,  R (T)  is  determined  for  two  temperatures,  Tt  and 
TH  say,  and  then  equation  3.7  and  its  counterpart  for 
perpendicular  components 


4 

Wf 


(2) 


(3.H) 


i 


can  be  evaluated  at  Tt  and  T„  and  combined  with  equation 

i_»  H 
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3-10  to  yield 


r(2)  . x(2)/x(2) 


He(TL)Re(TH)(TH-TL) 

TLTH[Re(TH)-Re(TL)]/Q+[THRe(TH)-TLRe(TL)] 


t(3)  . x(3)/x(3) 


(3.12) 


Re(TL)Re(TH)(TH-TL) 

Q[Re(TL)-Re(TH)]+[THR6(TL)-TLRe(TH)] 


where  Q,  defined  by 


(3.13) 


(3.14) 


is  minus  the  inverse  of  the  1/T  intercept  of  the  best-fit 

line  to  x6  vs.  1/T,  and  has  thus  been  determined  by  the 
II 

experiments  described  above.  Both  of  these  methods  were 
used  on  one  gas  as  a test  of  consistency,  and  the  results 
were  found  to  agree  within  experimental  uncertainty. 

3. 2d  |Ako| 

As  mentioned  in  section  3.1b  |Ako|,the  magnitude  of 
the  wave  vector  mismatch  defined  in  equation  2.29  must  be 
known  in  order  to  be  able  to  derive  x6.  |AkQ|  is  determined 
experimentally  by  observing  interference  between  second- 
harmonic  waves  which  are  generated  in  two  crystalline  quartz 


plates  by  a single  collimated  laser  beam,  in  an  experiment 
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devised  by  Finn  and  Ward.  A cell  placed  between  the  two 
plates  contains  the  vapor  under  investigation.  No  dc  field 
is  applied,  and  the  beam  is  not  focused,  so  no  second 
harmonic  is  generated  in  the  gas.  The  only  harmonic  fields 
present  are  those  generated  in  the  quartz  plates,  and  these 
fields  do  not  change  in  amplitude  with  changing  gas  density. 
However,  dispersion  in  the  gas  changes  the  relative  phase 
of  the  harmonic  generated  at  the  two  plates  as  gas  density 
changes.  Since  the  second  harmonic  leaving  the  second  plate 
is  the  sum  of  the  harmonic  waves  generated  at  the  two  plates 
its  intensity  varies  periodically  as  a function  of  the 
density.  It  is  assumed  here  that  n^  and  n^,  the  indices 
of  refraction  at  the  respective  frequencies,  are  linear 
functions  of  density,  which  is  a good  approximation  for 
dilute  gases. ^ The  magnitude  of  the  wave  vector  mismatch 
is  then  given  by^® 


| Ak  | = 
' o' 


2ttM 

DAp 


(3.15) 


where  Ap  is  the  density  difference  between  two  accurately 
located  harmonic  minima,  M is  the  number  of  intensity 
maxima  between  the  two  minima,  and  D is  the  distance  of 
experimental  gas  which  the  beam  traverses  between  the  plates 
Many  of  the  vapors  investigated  have  departures  from 
ideality  of  > 5#  at  the  temperatures  and  pressures  employed, 
so  it  is  important  to  evaluate  the  pressure  dependence  of 
the  density  correctly.  We  have  relied  upon  published  virial 
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coefficients  in  performing  such  evaluations. 

3. 3 Data  Analysis  and  Typical  Data 

The  analysis  of  |x6|  data  is  best  understood  by 
examining  a typical  data  run.  One  such  run  is  illustrated 
in  figure  5.  The  top  part  of  the  figure  is  a plot  of  the 
harmonic  versus  monitor  signals  for  each  laser  shot  of  a 
single  subrun.  The  dotted  lines  indicate  the  baseline 
setting  for  the  respective  signals.  The  value  of  the 
signal  S,  which  is  the  expected  value  of  the  ratio  of  the 
harmonic  to  monitor  channel  signals,  is  given  by 

N N 

S = E (H.-B  ) / E (M.-B  ) (3.16) 

i=l  1 H j=l  <:1  M 

where  is  the  number  of  counts  in  the  harmonic  channel 

for  shot  i,  Mj  is  the  harmonic  channel  count  for  shot  j, 

BH  and  BM  are  respectively  the  harmonic  and  monitor 

channel  baseline  settings,  and  N is  the  number  of  laser 

shots  in  the  subrun.  Equation  3.16  is  the  best  linear 

least-squares  fit  to  the  data  passing  through  the  point 

(Bm,Bh)  providing  that  all  of  the  statistical  fluctuation 

is  due  to  photon  (or,  more  accurately,  photoelectron) 

2 

statistics  in  the  monitor  channel,  i.e.,  that  cr^  is 


proportional  to  M.  (since  signal  levels  are  adjusted  so 
that  *■  H^,  assignment  of  fluctuations  to  the  monitor 
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(3.17) 


Ideally,  if  one  knew  the  gains  G^  and  of  the  harmonic 
and  monitor  channels,  such  that  the  number  of  photo- 
electrons emerging  from  the  harmonic  PMT  photocathode  was 
(Hi-B^/GR,  and  similarly  for  GM,  one  could  write 


“ N 

E (H.-B  ) 
i=l  1 M 


G, 


M 


(3.18) 


which  is  the  correct  variance  for  equation  3-16.  Such  a 
calibration  was  made  for  several  subruns,  and  it  was 
found  that  equations  3.17  and  3.18  were  in  good  agreement, 
indicating  that  the  data  spread  over  one  subrun  is  adequately 
accounted  for  by  photoelectron  statistics.  It  should  be 
noted  that  equations  3.16-3.18  are  based  upon  a Gaussian 
distribution  of  pulse  energies  for  a given  number  of  photo- 
electrons. The  |xe|  experiment  has  a conversion  efficiency 
of  approximately  10-12  to  10-1^  and  losses,  including  PMT 

quantum  efficiency,  of  around  a factor  of  100,  so  that  with 
17 

10  incident  laser  photons,  about  10  to  1000  photoelectrons 
leave  the  photocathode.  This  is  a sufficient  number  to 
satisfy  Gaussian  statistics,  at  least  approximately,  and  to 
be  insensitive  to  the  details  of  the  single-photoelectron 
output  distribution  of  the  photomultiplier. 
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The  bottom  part  of  figure  5 is  a plot  of  a typical 
|X6|  experimental  run.  Each  point  is  the  average,  in  the 
sense  of  equation  3.16,  of  a subrun  similar  to  that  in  the 
top  of  figure  5.  The  error  flags  are  the  n of  equation 

o 

3.17.  The  horizontal  scale  is  the  approximate  time  scale 

of  the  data  run.  Examination  of  this  run  reveals  what 

appears  to  be  a drift  of  the  signal  with  time,  outside  the 

limits  of  statistical  fluctuation.  Most  |xe|  runs  exhibit 

similar  drifts,  and  they  appear  to  be  due  to  slow  changes 

in  the  transverse  mode  structure  of  the  laser  beam,  drift 

of  the  beam  between  areas  of  different  quantum  efficiency 

on  the  photocathode  of  the  detector,  and  drift  of  the  beam 

focus  in  the  region  of  the  dc  field.  The  beam  is  known10 

to  have  components  of  Gaussian  modes  of  other  than  the 

lowest  order,  and  the  dcSHG  process  is  quite  sensitive  to 

iq 

the  transverse  mode  structure.  Alternating  subruns  of 

experimental  gas  with  those  of  reference  gas  has  the  effect 

of  recalibrating  C in  equation  3-3  at  brief  intervals, 

minimizing  the  effects  of  the  drifts,  and  run-to-run 

fluctuations  are  consistent  with  the  subrun  statistics  of 

equation  3-17.  An  absolute  calibration  of  C,  including 

analysis  of  transverse  mode  structure,  evaluation  of 

losses  in  optical  components  and  gains  and  losses  in 

detectors  and  electronics,  and  measurement  of  the 

efficiency  of  the  quartz  harmonic  generator  in  the  monitor 

5 6 

box,  yields  absolute  measurements  of  x with  uncertainties 
of  about  a factor  of  3*  This  large  uncertainty  precludes 
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useful  comparison  with  theory.  Relative  measurements  of 
0 

X for  pairs  of  molecules,  on  the  other  hand,  yield  results 
with  much  smaller  uncertainties  in  the  range  of  1-5$. 

We  put  these  results  on  an  absolute  basis  without  intro- 
ducing a large  experimental  uncertainty  by  using  an  accurate 
calculated  value  for  helium  as  described  in  section  3.2a. 

Thus,  the  information  obtained  from  the  run  illus- 
trated in  figure  5 is  S(T)  for  the  experimental  gas  and 
S’  for  the  reference  gas,  which  are  the  averages  of  the 
respective  quantities  from  the  individual  subruns.  Voltages 
are  adjusted  so  that  S(T)  and  S'  are  nearly  equal,  minimizing 
systematic  errors  due  to  baseline  inaccuracies.  Typical 
uncertainty  of  the  ratio  of  |xe/AkQ|  of  the  experimental 
gas  to  that  of  the  reference  gas  is  1$. 

Figure  6 shows  a typical  plot  of  x^  against  inverse 
temperature.  The  error  bars  are  comparable  with  the  size 
of  the  data  points,  and  include  the  uncertainty  of  |AkQ| 
(where  applicable,  the  uncertainty  of  x6  would  also  include 

( ■7) 

that  of  Xj | of  methane  relative  to  that  of  helium,  but 

the  uncertainty  of  the  calculated  value  for  helium  is  not 

included).  The  line  through  the  points  is  the  best  least- 

C (o) 

squares  fit  to  the  data,  and  xv  ' and  xv  are  derived 

il  ti 

using  equations  3.8  and  3*9>  respectively,  with  p.  obtained 

(2) 

from  the  literature.  Uncertainties  are  determined  for  xv 

( *5  \ e 

and  xv  from  the  uncertainties  of  the  x and  4 data  by 

20 

standard  methods. 


0 

R (T)  is  determined  with  a statistical  accuracy  of 

1 to  3#.  Systematic  effects  in  Re(T)  are  discussed  in 

reference  17,  and  are  small.  In  the  cases  where  the 

slope  of  x6  vs.  1/T  is  small,  there  is  a large  relative 

(2) 

error  in  Q,  and  Rv  ' is  ill-determined.  In  these  cases, 

(2) 

Rv  ’ is  not  reported,  and  R'J'  and  its  uncertainty  are 

(2) 

calculated  with  the  assumption  that  Rv  ' is  equal  to 
3.0+0. 3,  consistent  with  the  better  determined  values. 

Uncertainties  in  AkQ  are  typically  1$.  This  includes 
an  approximately  0.5#  uncertainty  in  the  distance  Ap  between 
measured  intensity  minima  (a  combination  of  statistical 
and  pressure  measurement  uncertainties)  and  5#-10# 
assumed  uncertainties  in  the  second  virial  coefficients 
of  the  gases  as  determined  from  the  literature. 

We  have  seen  how  the  dcSHG  theory  described  in 
chapter  2 is  implemented  in  measuring  effective  third-order 
polarizabilities  in  gases.  Further,  we  have  seen  how  the 
theory,  when  applied  to  such  measurements  made  as  a function 
of  temperature  and  polarization,  along  with  dispersion 
measurements  and  other  molecular  data  can  be  used  to 
determine  second-  and  third-order  polarizabilities  of 
molecules.  In  the  next  two  chapters,  we  will  see  the 
results  of  these  measurements  for  several  molecular  gases. 
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CHAPTER  4 

EXPERIMENTAL  RESULTS  FOR  ASSORTED  SMALL  MOLECULES 

We  have  applied  the  experimental  techniques  presented 
in  chapter  3 to  a number  of  molecules.  These  can  be 
divided  into  two  broad  categories  — chloro-  and  fluoro- 
substituted  (halogenated)  methanes,  and  "other"  molecules. 
The  halogenated  methane  data  will  be  presented  in  chapter 
5,  and  will  be  further  examined  in  conjunction  with  the 
model  discussed  in  chapter  6. 

The  assorted  "other"  molecules  were  selected  for  a 
number  of  reasons;  data  from  other  nonlinear  optical 
experiments  exist  for  several;  theoretical  results  are 
available  for  a few;  and  interesting  correlations  can  be 
applied  to  some  groups  of  them.  All  of  these  molecules 
share  the  property  of  being  relatively  small,  simple 
systems  for  which  hyperpolarizability  calculations  in  the 
domain  of  molecular  orbital  theory  are  feasible. 

The  experimental  parameters  and  data  for  these  gases 
are  presented  in  section  1.  Also  shown  there  are  the 
derived  hyperpolarizabilities  and  ratios,  along  with  other 
experimental  hyperpolarizability  data  and  theoretical 
estimates.  In  section  2,  special  points  of  interest  for 
certain  of  the  molecules  are  discussed. 
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4. 1 Experimental  Parameters,  Data,  and  Derived  Quantities 

The  experimental  parameters  and  data  for  the  molecules 

not  possessing  permanent  dipole  moments,  and  therefore 
0 

having  x which  are  independent  of  temperature,  are  dis- 
played in  table  1.  The  values  of  E°  are  the  field  strengths 
on  the  z axis  between  the  electrodes,  as  calculated  from 
equation  2.26  with  z = 0,  for  typical  voltages  V . Values 
of  AkQ  were  all  negative  for  these  gases,  and  were  obtained 
from  experimental  data  with  the  aid  of  the  virial  data 

referenced  in  the  table,  x6  and  Re  were  measured  at  room 

II 

temperature  (approximately  20°C)  and  optimum  density  0OT)t- 

The  uncertainties  listed  for  X^  are  standard  deviations. 

Similar  parameters  and  data  for  the  dipolar  molecules 

0 

are  shown  in  table  2.  The  temperature  Tq  at  which  R was 
measured  is  also  listed  for  these  molecules.  In  addition, 
since  x*  and  Re  are  both  temperature-dependent,  their 
slopes  with  respect  to  1/T  and  their  inf inite-temperature 
limits  are  indicated.  The  first  method  of  section  3.2c 
was  used  to  determine  R^2^  and  R^^,  so  that  only  one 
value  of  Re  appears  for  each  gas.  In  the  cases  of  CO,  NO, 
and  HpS,  there  was  Insufficient  signal  to  recover  statistic- 
ally meaningful  values  for  x^;  however,  for  these  gases  Re 
did  not  differ  significantly  from  3.0  at  room  temperature, 
and  there  was  no  indication  of  unusual  behavior  of  x^. 

The  values  of  xjj^  and  R^^  for  both  dipolar  and 
non-dipolar  molecules  as  derived  from  the  experimental 
data  are  shown  in  table  3,  along  with  results  from  other 
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Table  1 . 


E° 

( esu) 

°opt 

-“oa 

(-cm"1) 

1036xe 

II 

(esu) 

Re 

H2 

35 

2.97 

1.27+.01 

.06521.0008 

2. 861. 03 

n2 

36 

1.86 

1.89+. 02 

.08661.0010 

3.001.06 

°2 

48 

1.46 

2.4l±.02 

.09531.0016 

2.971.07 

o 

o 

ro 

43 

1.05 

3. 361. 03 

. 11191.0013 

2. 811.06 

SP6 

58 

1. 08 

3.471.03 

.130  1.002 

2.971.07 

a)  Using  virial  data  from: 

Reference  42  (H3,  Na , Oa ) ; K.  E.  MacCormack  and  W.  G. 
Schneider,  J.  Chem.  Phys.  18,  1269  (1950)  (C0a);  and 
H.  P.  Clegg,  J.  S.  Rowlinson,  and  J.  R.  Sutton,  Trans. 
Faraday  Soc.  51,  1327  (1955)  (SFs). 
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Table  3. 


■gj 

■E 

Kerr 

THGb 

Theory 

HP 

2.86 
± .03 

.0652 

±.0008 

• 079d 

.047® 

±.005 

.080 

±.012 

.034h 

N2 

3.00 
± .06 

. 0866 
±.0010 

.104 

.12  f 
±.01 

.107 

±.017 

. 0711 

°2 

3.12 
± .07 

.0953 

i.0016 

.11  d 

. 100 

sf6 

2.97 
± .07 

.130 

±.002 

.20 

.192 

.20  g 
±.02 

co2 

2.81 
± .08 

.1119 

±.0013 

.192 

.75  f 
±.16 

.156 

±.023 

CO 

2.92C 
* .15 

.144 

±.004 

.138 

NO 

3.32° 

± .15 

.235 

±.007 

.322 

h9s 

2.84c 
± .11 

.865 

1.022 

nh3 

3.19 
± .11 

+•  511 
±.009 

h2o 

.2.94 

- .55 

.194 

±.010 

(ch3)2o 

2.99 
± .15 

.529 

±.011 

CH30H 

2.90 
± .21 

.385 

±.011 

a)  Three-wave  mixing  (TWM)  data  from  reference  43 
normalized  using  the  dcSHG  coefficient  for  argon  from 
reference  13. 

b)  Third-harmonic  generation  ( THG)  data  from  J.F.  Ward  and 
G.H.C.  New,  Phys.  Rev.  185,  57  (1969). 

c)  Derived  assuming  Rv3)  JT0±0.3. 

d)  Data  from  reference  4,  normalized  using  the  dcSHG 
coefficient  for  argon  from  reference  13. 

e)  A. D.  Buckingham  and  B.J.  Orr,  Proc.  Roy.  Soc.  A 305, 

259  (1968). 


f)  Reference  21. 

g)  A.D. 


g)  A.D.  Buckingham  and  D. A.  Dunmur,  Trans.  Faraday  Soc. 
64,  1776  (1968). 

h)  Xzzzz (° ;0#0,0)  from  reference  44. 


i)  X 


zzzz 


(0;0,0,0)  from  reference  22. 
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gas-phase  experiments,  and  available  theoretical  estimates. 

does  not  differ  significantly  from  3-0*  the  value 
which  would  apply  in  the  limit  of  small  dispersion  in  the 
molecular  This  result  is  not  surprising  since  all  of 

the  experimental  frequencies  involved  are  far  from  molecular 

resonances.  The  other  experimental  gas-phase  data  for 

fq') 

Xjj  are  from  other  dcSHG  experiments  and  from  three-wave 

mixing  ( TWM ) , Kerr  effect,  and  third-harmonic  generation 

( THG)  experiments.  These  data  tend  to  be  slightly  larger 

than  our  values,  but  not  unreasonably  so  in  view  of 

experimental  uncertainties  and  the  likelihood  of  small 

effects  due  to  dispersion.  A possible  exception  is  x^^ 

for  CO^j  for  which  the  Kerr  effect  measurement  of 

21 

Buckingham  et  al.  is  much  larger  than  other  values, 
whereas  their  data  for  N2  and  SFg  agree  acceptably  with  the 
others. 

The  theoretical  values  of  xjj^  listed  for  and  N2 
are  the  molecular  \ only.  The  value  listed  for  N0  is 

ZZZZ  2 

22 

preliminary,  and  it  is  anticipated  that  the  calculation 
will  be  refined  and  extended  to  the  other  components  of 
*i1kl  for  tlle  8ases  investigated  in  this  dissertation  and 
in  reference  12. 

(2)  (2) 

Values  for  xv  ' and  Rv  ',  as  derived  from  the 
II 

experimental  data  and  from  literature  values  for  permanent 
dipole  moments,  are  shown  in  table  4,  along  with  other 
experimental  data  and  with  theoretical  results.  Comments 
on  the  single  other  experiment,  a dcSHG  measurement  in 
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Footnotes  for  Table  4. 

a)  Electric  dipole  moment  in  Debye  units  from  references 
45  and  46. 

b)  N.  S.  Hush  and  M.  L.  Williams,  Theoret.  Chim.  Acta 
(Berlin)  25,  346  (1972). 

c)  x^2z  (0;0,0)  from  reference  44. 

d)  Reference  23. 

e)  S.  P.  Liebmann  and  J.  W.  Moskowitz,  J.  Chem.  Phys.  54, 

3622  (1971).  — 

f)  Reference  24. 

g)  Reference  25. 

h)  Refractivity  virial  data  from  reference  47. 
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liquid  H^O,  is  deferred  to  the  next  section. 

The  theoretical  results  are  from  molecular  orbital 

calculations,  and  are  listed  according  to  the  general 

method  employed,  either  "semi-empirical"  (SE),  uncoupled 

Hartree-Fock  (HF),  or  coupled  Hartree-Fock  (CHF).  Wide 

discrepancies  in  magnitude,  and  even  sign  differences, 

occur  between  the  results  of  different  authors  using 

various  techniques.  For  example,  in  H^O,  the  theoretical 
•?0  (2) 

results  for  10  1 range  from  -0.12  to  +0.09  and  in  CO, 

they  vary  from  -0.44  to  +0.88.  In  addition,  results  by 

the  same  authors  vary  considerably  depending  on  the 

2^  24  2E 

number  and  nature  of  the  basis  functions  used.  * 

Our  data,  with  an  uncertainty  of  less  than  21$  in  the  worst 
case,  and  3-12$  for  the  rest,  can  provide  a stringent  test 
of'  molecular  orbital  calculations  for  these  molecules. 

4. 2 Discussion  of  Results  for  Particular  Molecules 

Certain  of  the  12  molecules  considered  in  this  chapter 
warrant  particular  attention,  and  they  are  discussed  in 
this  section  in  turn. 

4.2a  H2,  N2,  and  02 

Because  of  the  relative  simplicity  of  the  homonuclear 
diatomic  molecule,  it  seems  reasonable  to  attempt  to 
predict  X(3)  for  H2,  N2,  and  02  using  an  adaptation  of  an 

p/T 

approximation  due  to  Dawes  which  meets  with  some  success 
in  predicting  the  variation  of  among  the  inert 

gases.  -*-3*26  In  pawes  approximation,  the  perturbation 
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f 


expansion  for  obtained  by  summing  terms  such  as 

that  in  equation  2.5,  is  reduced  to  a single  term  in  a 
manner  analogous  to  the  linear  polarizability.  A single 
effective  state  frequency  a>Q  is  substituted  for  transition 


freauencies  ax,  . in  order  to  factor  out  the  term  a> 
i j o 

the  approximation  is  made  that 
-1 


-5 


Next , 


/ 2V 

< r > = 03 

gg  o 


E co  . ( r) 


( r > 


gi  - gi  - ig 


(4.1) 


where  g refers  to  the  ground  state.  Finally,  the 

assumption  is  made  that  is  the  sum  of  Nq  identical 

"one  electron"  expressions,  where  Nq  is  an  effective 

electron  number.  The  resulting  expression  is  proportional 

to  a sum  like  that  on  the  right-hand  side  of  equation  4.1, 

27 

and  the  Thomas-Reiche-Kuhn  sum  rule  is  invoked  to  collapse 

it  to  a single  term.  The  N and  co  are  then  taken  to  be 

o o 

the  values  used  to  fit  the  dispersion  of  the  linear  index 
of  refraction.  In  the  case  of  homopolar  diatomic  molecules, 
one  can  use  data  on  the  anistropy  of  the  linear  molecular 
polarizability  to  fit  two  effective  electron  numbers  (which 

can  be  thought  of  as  the  numbers  of  electrons  oscillating 
respectively  parallel  to  and  perpendicular  to  the  inter- 
nuclear  axis).  The  expression  arrived  at  after  averaging 
over  molecular  orientations  is 
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Li 


1 


0 

0 

D 


N e 

x(3)  = — 1L — _ 
" hm2<n5 


(l-4r2)2(l-r2)2 


(4.2) 


in  which  the  quantities  are  as  follows:  6 is  1-(N  /N  ), 

where  N and  N are  the  effective  electron  numbers  in  the 

ll  a 

parallel  and  perpendicular  directions,  respectively;  r is 

the  ratio  's>/'x>Q,  where  cd  is  the  field  frequency;  and  e is 

the  charge,  and  m the  mass,  of  the  electron.  Since  the 

Dawes  approximation  reproduces  the  values  of  for  the 

inert  gases  only  to  an  overall  multiplicative  factor,  it  is 

necessary  to  normalize  the  Dawes  values  of  x ^ ^ for  the 

II 

homopolar  diatomic  molecules  to  that  for  helium.  Table  5 

shows  the  results  of  the  application  of  equation  4.2  to 

H0,  Nn,  and  0o,  with  values  of  N and  cd  obtained  from 
d d d o o 

reference  28  and  polarizability  anisotropy  data  from 

reference  29.  While  far  from  a good  fit,  the  Dawes 

approximation  values  are  of  the  correct  order  of  magnitude, 

and  display  the  correct  direction  of  variation  from 

( q \ 

molecule  to  molecule.  Even  when  normalized  to  x of  N^, 
the  Dawes  values  for  and  0^  do  not  agree  convincingly 
with  the  experimental  values.  However,  because  of  the 
severity  of  the  approximations,  and  the  increased  complexity 
of  homonuclear  diatomic  molecules  over  that  of  inert  gas 
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atoms,  it  is  not  astonishing  that  the  Dawes  approximation 
should  be  less  applicable  to  these  molecules. 

4.2b  NH3 

The  permanent  electric  dipole  moment  of  the  static 
pyramidal  NH^  configuration  was  used  in  table  4.  The 
ground-vibrational  state  of  NH^  is  inversion-split  by 
the  energy  e corresponding  to  the  frequency  of  tunneling 
of  the  N atom  through  the  center  of  the  triangle  formed  by 
the  three  H nuclei.  This  splitting  is  small,  however,  and 
it  can  be  shown  that  for 

e2/(kT)2  « 1 (4.3) 

ammonia  behaves  like  a gas  of  ordinary  dipolar  molecules. 

2 2 -S 

At  room  temperature  the  value  of  e /(kT)  is  about  2x10  , 

so  the  inequality  4.3  is  well  satisfied.  Indeed,  the  fact 

that  x is  a linear  function  of  1/T  in  the  range  of 

temperatures  used  here  (see  figure  5)  confirms  the  ordinary 

behavior  of  ammonia. 

The  ammonia  molecule  also  illustrates  dramatically  the 

fact  that  the  differences  of  R^2^  and  R^)  from  3*0,  while 

small,  are  indeed  non-zero.  If  R^2^  and  r(3)  were  equal, 

then  R (T)  would  have  to  be  a constant  function  of 

temperature,  equal  to  R^2^  (and  to  R^'),  as  can  be  seen 

from  equations  3.7  and  3.10-3.13.  However,  as  shown  in 

tables  2,  3,  and  4,  Re  is  significantly  different  from 
(2)  (R) 

Rv  ' and  Rw/,  so  that  the  latter  two  quantities  cannot 
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be  equal  (in  particular,  they  cannot  both  be  equal  to  3.0) 
indicating  that  the  effects  of  dispersion  in  various 
components  of  Xj_jk  and  while  they  might  be  small, 

cannot  be  neglected  entirely. 

4.2c  H20 

•50  e 

Levine  and  Bethea  have  measured  x in  water  and  in 

II 

water-methanol  and  water-dioxane  solutions  by  the  method 
of  dcSHG  in  the  liquid  phase.  Because  the  details  of 
intermolecular  interactions  in  strongly  associating  liquids 
are  not  well  understood,^  Levine  and  Bethea  were  not  able 

to  extract  precise  isolated-molecule  values  of  and 

li 

( -3  1 

Xjj  ' for  water  from  their  data.  By  postulating  a model 

of  the  interactions,  however,  they  arrive  at  the  expression 
(2) 

for  x , the  second-order  polarizability  per  molecule  of 

s 


water  in  solutions,  of 


(2)  _ 

■s 


xr;  = (XT'  + x 


(2)  + v(2) 


hb  ) + xdip(1_g) 


(4.4) 


(2) 

where  xv  ' is  the  value  for  an  isolated  molecule 
II 

/ p \ 

(corresponding  to  our  gas-phase  value),  X^^  is  the 
second  order  nonlinearity  due  to  hydrogen  bonding  among 

is  the  effective 


the  water  and  solvent  molecules,  x 


dip 


second-order  hyperpolarizability  induced  in  the  molecule 
by  the  interaction  of  the  molecular  x^  with  the  average 
electric  field  due  to  the  permanent  dipole  moments  of 
surrounding  molecules,  and  g is  the  Kirkwood-Frolich 
correlation  parameter,  having  to  do  with  short-range 


» i 
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( O \ 

molecular  alignment.  They  derive  values  of  xv  from 

s 

their  data  using  an  estimate  for  the  gas-phase  value  of 

Xv  ' consistent  with  our  measured  value.  They  calculate 
II 

values  of  g for  mixtures  from  data  on  the  dielectric 

(2) 

constants  of  mixtures.  Their  derived  values  of  x v 

^ s 

range  from  about  1/2  to  1/8  of  our  x^2^  value,  and  they 

(2)  11 

find  that  in  this  range,  x~  ' varies  approximately  as 

o 

(1-g).  This  behavior  along  with  equation  4.4,  implies  that 


\/)+*hb*°  <4-5) 

to  within  the  broad  uncertainty  limits  of  the  data  and 
model.  Such  a cancellation  would  be  interesting,  in  that 
it  would  suggest  that  hydrogen  bonding  might  cause  H and 
0 atoms  to  be  arranged  with  a high  degree  of  short-range 
symmetry.  For  example,  an  0 atom  might  be  surrounded  by 
four  H atoms  in  a T^  symmetry,  while  an  H atom  could  lie 
between  two  0 atoms  in  a C symmetry.  The  value  of 

® V 

(2) 

Xv  ' , as  defined  by  equation  2.21,  is  zero  for  both  of 

p 

these  symmetries. 


In  this  chapter  we  have  presented  experimental 

parameters  and  data  for  \e  and  Re,  and  their  temperature 

dependence,  from  which  we  have  extracted  R^, 

(2)  ii  li 

and  Rv  ' values  for  a number  of  simple  molecules.  Our 

data  are  largely  consistent  with  other  experimental  values 

that  are  available  for  comparison.  Where  the  values 

disagree,  we  feel  that  our  data  are  less  subject  to  errors 


A 
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in  measurement  and  interpretation.  Our  experiments  have 
an  inherently  high  signal-to-noise  ratio,  and  systematic 
errors  are  compensated  for  in  a straightforward  manner, 
with  calibration  an  integral  part  of  the  experiment.  In 
addition,  the  molecular  observables  are  derived  from  the 
experimental  data  in  an  unambiguous  manner,  without 
having  to  estimate  contributions  to  the  effects  either 
from  other  processes  or  from  intermolecular  interactions. 
While  we  have  commented  upon  particular  features  of  some 


of  the  molecules  here,  we  defer  a discussion  of  water, 
methanol,  and  dimethyl  ether  as  a group  to  the  section 
on  the  Bond  Additivity  Approximation  in  chapter  6. 


CHAPTER  5 

EXPERIMENTAL  RESULTS  FOR  SOME  HALOGENATED  METHANES 


1 

I 


The  halogenated  methanes  investigated  here  are  all 

of  the  form  CX  , with  n = 1-4  and  X,Y  = F,C1,H.  The 

(2\  (2) 

second-  and  third-order  polarizabilities  xv  ard  xv  for 

some  of  the  molecules  have  been  measured  by  Bigio  and 
12 

Ward.  The  capability  of  making  gas-phase  measurements  in 
substances  which  are  liquids  at  room  temperature  has 

1 2) 

allowed  us  to  complete  the  series  with  measurements  of  xv 

and  X(3)  for  CCl^,  CHCl^  CH2C12,  CHgCl,  noF2  and  CC^F. 

In  addition,  we  have  made  an  especially  ' se  measurement 

(1) 

of  xjj  ' for  CH^,  since  it  was  used  here  a erence  gas. 

The  data  presented  in  this  chapter  have  appeared  in  The 

Physical  Review,  Section  A.3^- 

The  experimental  parameters,  and  experimental  data 

on  dispersion,  effective  hyperpolarizabilities,  and  effective 

ratios,  are  shown  in  table  6.  These  data  are  similar  in 

nature  to  those  presented  in  tables  1 and  2.  However,  for 

these  molecules,  the  second  method  discussed  in  section  3.2c 

(2)  ( q ) 

for  determining  Rv  ' and  RVJ/  was  used.  Thus,  there  are 

two  values  listed  for  Re,  as  well  as  two  values  for  the 

temperatures  at  which  the  Re  were  measured,  for  each  gas 

(except  for  the  non-dipolar  molecules  CH^  and  CCl^,  for 

which  a single  value  of  Re  suffices). 

The  results  for  x^3^  and  are  shown  in  table  7, 

II 

along  with  other  gas-phase  experimental  data.  Again,  R^) 
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Table  7. 


r(3) 

y(3) 

*11 

Other  x(3)  data 
dcSHG  tIIm  Kerr 

CC1U 

3.02 

1.37 

1.7  C 

1.7  C 

1. 

66  e 

±0.09 

±0.02 

±0. 

05 

CHC1-3 

2.96 

1.12 

1.19c 

3 

±0.05 

±0.03 

0HoCl_ 

3>°2 

0.92 

2 2 

±0.09 

±0.03 

CH_C1 

3.28 

0.57 

3 

11.03 

±0.03 

ch4 

3.°5b 

±0.05 

0.263 

±0.003 

0.32C 

0.34= 

0.19 

0. 

±0. 

242  r 
0121 

CC10F0 

3.55 

0.61 

±0.35 

±0.01 

CCl-F 

3.39 

0.96 

5 

12.92 

*0.03 

al  Using  the  assumption  r(2)  = 3. 0-0-3. 

b)  From  reference  12. 

c)  Data  from  reference  4 normalized  using  the  dcSHG 
coefficient  for  argon  from  reference  13. 

d)  Data  from  reference  43  normalized  using  the  dcSHG 
coefficient  for  argon  from  reference  13. 

e)  Reference  48  . 

f)  A. D.  Buckingham  and  B.J.  Orr,  Trans.  Faraday  Soc.  65, 
673  (1969). 
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is  not  significantly  different  from  3.0,  and  as  before, 

(2) 

for  cases  in  which  Rv  ' cannot  be  accurately  determined, 

r(3)  calculated  by  assuming  the  value  R^2 * * *^  = 3.O+O.3. 

Other  experimental  data  are  from  TWM,  THG,  and  Kerr  effect 

measurements,  and  there  is  excellent  agreement  between 

these  data  and  ours.  A number  of  measurements  have  been 

made  of  the  third-order  nonlinearity  of  CCl^  in  the  liquid 

phase  for  various  processes,  such  as  dcSHG,  Kerr  effect, 

and  ellipse  rotation.  Much  of  this  information  has  been 

32 

reviewed  by  Levine  and  Bethea,  who  use  Onsager  local- 

q ( q\ 

field  factors^  to  extract  estimates  for  x j ( of  the 

isolated  molecule.  They  arrive  at  estimates  which  are 

typically  within  15%  of  our  experimental  result.  There 

are  currently  no  at)  initio  calculations  of  available 

for  these  molecules. 

(2)  (2) 

Results  for  x'  and  Rv  are  shown  in  table  8,  along 

(2) 

with  dipole  moments  and  measurements  of  x v from  other 

(2)  11 

sources.  Rv  ' is  subject  to  large  uncertainties,  and  in 

fact  could  not  be  assigned  a meaningful  value  in  the  case 

of  CHCl^  and  CH^Cl^  (for  which  the  change  in  x6  over 

the  temperature  range  used  was  very  small).  For  the  other 
(2) 

molecules,  Rv  ' does  not  differ  significantly  from  3.0. 

(2) 

The  sign  of  x|  > assigned  according  to  the  conventions  of 

section  3.8b,  is  negative  for  CC12F2  and  CCl^F,  as  it  is12 

for  CCIF^  and  for  CHnF^  n,  n = 1-3,  but  it  is  positive  for 

CHnCl^_n  (however,  negative  and  zero  values  of  y/2'' 
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within  the  uncertainty  limits  for  CHCl^).  The  uncertainties 

of  other  published  gas-phase  measurements  of  xj2^  are 

large,  so  although  the  other  data  tend  to  be  larger  in 

magnitude  than  our  values,  they  are  not  necessarily 

inconsistent.  A possible  exception  is  the  value  from  the 

4 

dcSHG  and  TWM  experiments  of  Mayer  £t  al.  for  CHCl^,  which 

differs  from  our  result  by  40  times  our  experimental 

(o\ 

uncertainty.  It  should  be  noted  that  xjj  ' is  obtained  by 
Mayer  e_t  al.  from  the  difference  between  the  results  of  the 
dcSHG  and  TWM  experiments,  which  require  independent 
calibrations  with  different  sets  of  frequencies,  whereas 
in  our  method,  only  relative  measurements  of  signal  as 
a function  of  temperature  are  made,  all  at  the  same  set  of 

qo  (2) 

frequencies.  If  10J  Xjj  were  10  times  our  uncertainty, 
or  approximately  0.1  esu,  there  would  have  been  an  easily- 
seen  15$  change  in  our  experimental  signal  over  the 
temperature  range  319-500°K,  rather  than  the  ~ 1$  change 
actually  observed. 

The  methods  for  extracting  isolated-molecule  values  of 

(2) 

Xv  from  measurements  in  the  liquid  phase  are  indirect, 
involving  estimation  of  poorly-understood  effects  due  to 
intermolecular  interactions.  The  agreement  between  these 
data  and  our  measurements  is  poor. 

(2) 

There  are  currently  no  ab  initio  calculations  of  X 
for  any  of  these  halogenated  methanes. 


We  have  presented  experimentally  determined  temperature- 
dependent  values  of  x^  and  Re  for  some  halogenated  methanes, 
as  well  as  the  derived  X^,  R^3^,  and  R^2)  values, 

and  compared  them  with  other  available  experimental  data. 
Again,  our  values  are  for  the  most  part  consistent  with 
other  data,  and  where  this  is  not  the  case,  we  feel  that 
our  values  are  the  more  reliable,  as  discussed  in  chapter  4. 
In  the  next  chapter,  we  will  investigate  correlation  of 
some  of  these  properties,  along  with  other  electric  tensor 
properties,  among  the  molecules  of  the  complete  halogenated 
methane  series. 
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CHAPTER  6 

CORRELATION  OF  ELECTRIC  TENSOR  PROPERTIES 
AMONG  SIMILAR  MOLECULES 

(2) 

Molecular  electric  tensor  properties  a,  xv  ' , 

X<3),  and  so  on  are  ultimately  determined  by  the  changes 
in  electronic  wave  functions  due  to  applied  electric 
fields,  and  can  be  written  formally  in  terms 

7 

of  time-dependent  perturbation  theory  expressions.  Such 
expressions,  while  correct  in  principle, generally  defy 
evaluation  because  of  the  difficulty  of  finding  molecular 
wave  functions.  However,  the  obvious  structural  similarity 
among  the  molecules  within  certain  groups  suggests  that 
the  electric  tensor  properties  of  the  molecules  in  a 
given  group  might  be  correlated.  In  such  correlations, 
the  electric  tensor  properties  of  the  set  of  molecules 
are  calculated  from  a smaller  set  of  parameters.  For 
the  cases  discussed  here,  the  parameters  are  interpreted 
physically  as  the  electric  tensor  properties  of  the  bonds 
or  structural  groups  in  the  molecules,  and  the  predicted 
tensor  properties  for  a given  molecule  are  expressed  as  a 
function  of  the  properties  of  the  component  segments  and 
the  geometry  of  the  molecule. 

In  section  1,  we  examine  the  Bond  Additivity 
Approximation  (BAA),  in  which  the  electric  tensor 
properties  of  a molecule  are  assumed  to  be  the  sums  of 
respective  properties  associated  with  its  component  bonds 
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68 

or  structural  groups.  In  section  2,  an  Interacting 
Segment  Model  (ISM)  of  these  properties  is  proposed, 
with  the  aim  of  avoiding  some  of  the  drawbacks  of  the 
BAA  while  retaining  much  of  its  simplicity.  In  section 

3,  the  ISM  is  applied  to  the  permanent  dipole  moments 

4,  linear  polarizabilities  a,  and  second  order  polar- 

(2) 

izabilities  xv  of  the  halogenated  methanes  in  order  to 
obtain  a clearer  picture  of  the  contributions  of  bond 
interactions  to  the  electric  tensor  properties  of  this 
set  of  similar  molecules.  Much  of  the  material  in  this 

33 

chapter  has  appeared  in  The  Journal  of  Chemical  Physics. 

6 . 1 The  Bond  Additivity  Approximation 

The  concept  of  assigning  separate  electric  tensor 
properties  to  the  bonds  (such  as  the  C-Cl  bond  in  CH^Cl') 
and/or  functional  groups  (such  as  the  CH^j  or  methyl, 
group  in  CH^OH)  of  a molecule,  and  summing  the  respective 
properties  over  the  segments  (bonds  and  groups)  to 
determine  the  overall  properties  of  the  molecule,  is 
called  the  Bond  Additivity  Approximation.  This  model  has 
been  widely  (although  not  universally)  accepted  by 

•3h 

chemists  and  physicists  for  more  than  fifty  years, 
largely  because  of  its  simplicity.  The  BAA  has  found  use 
in,  for  example,  resolving  questions  about  molecular 
geometry. 

In  the  BAA,  the  ith  segment  in  a molecule  is 
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assigned  tensor  properties  x(n),  n = 0,1,..., 

is  the  permanent  electric  dipole  moment,  x^1) 
the  linear  polarizability,  and  so  on  (we  will  use 
and  x^  interchangeably  with  p.  and  a,  respectively,  as 
appropriate).  The  BAA  value  x/n^  for  the  molecular 
tensor  is  found  by  adding  the  respective  segment  tensors: 


(6.1) 


The  segment  tensors  are  here  assumed  to  be  independent  of 
their  environment  so  that,  for  example,  xjn^(C-Cl)  in 
CH^Cl  is  equal  to  x[n^(C-Cl)  in  CCl^  (in  the  coordinate 
system  of  the  bond).  In  actuality,  the  quantities  which 
are  to  be  reproduced  by  any  model  such  as  this  are  scalar 
observables  related  to  the  electric  tensors  (for  example, 

(2)  (2) 

we  know  xjj  , rather  than  x^^.).  The  relation  between 
each  tensor  property  and  its  associated  observable  will 
be  defined  when  it  is  introduced. 

We  have  two  sets  of  molecules  on  which  to  test  bond 
additivity.  The  first  set  has  only  two  segment  types, 
and  consists  of  only  three  molecules— water  (H20),  dimethyl 
ether  ((CH^)?0),  and  methanol  (CH^OH).  The  second  set 
has  three  types  of  segments,  and  consists  of  twelve 
halogenated  methanes.  We  will  discuss  the  smaller  set  of 
■ olecules  first  because  of  their  structural  simplicity. 


! 
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The  three  molecules  in  the  first  set  may  be  char- 
acterized as  an  oxygen  atom  with  two  groups  attached. 

The  groups  may  both  be  hydrogen  atoms,  or  both  CH^  groups, 
or  one  hydrogen  atom  and  one  methyl  group,  constituting 
water,  dimethyl  ether,  and  methanol,  respectively. 

We  will  write  the  properties  associated  with  the  G-H 
bond  as  xJn^(0-H),  and  those  associated  with  the  0-C  bond 
and  CH^  group  collectively  as  _x^n^  (0-CH^) . As  discussed 
later,  the  question  of  the  location  within  the  segment  at 
which  the  properties  act  is  not  important  in  the  BAA,  and 
will  be  ignored  here.  The  symmetry  axis  of  the  CH^  group 
is  considered  to  be  along  the  0-C  internuclear  axis,  so 
that  the  (O-CH^)  segment  symmetry  axis  will  also  be  in 
this  direction.  The  two  properties  of  these  molecules 
that  we  will  examine  first  are:  the  isotropic  linear 
polarizability  a,  defined  by 


a = -».(a  + a + a ) ; 

V xx  yy  z z ' * 

(6.2) 

and  which,  as  we  have  already  seen  in 

equation  2.19, 

is  a twice-contracted  fourth-rank  tensor. 

Both  of  these 

quantities  transform  as  scalars,  and  therefore  the  segment 

quantities  add  very  simply: 

a(H20)  = 2a(0-H) 

(6.3) 

a((CH3)20)  = 2a(0-CH^ 

(6.4) 

u 

(I 


I! 


f] 

n 
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a ( CH^OH)  = a(O-H)  + ct(0-CH3)  (6.5) 

which  means  that  BAA  predicts  the  isotropic  linear 
polarizability  per  molecule  of  methanol  to  be  the 
simple  average  of  those  of  water  and  dimethyl  ether. 

The  same  statement,  along  with  the  analogs  of  equations 


6. 3-6. 5,  apply  to  xjj3^: 

xj3)(H2°)  = 2XJ3)(0-H) 

(6.6) 

xJ3)((CH3)20)  = 2Xf|3)(0-CH3) 

(6.7) 

xJ3)(ch3oh)  = xft3)(o-H)  + x|(3)(o-ch3) 

(6.8) 

We  choose  a(O-H),  a(0-CH3),  x|(3^(0-H),  and  x 
to  satisfy  equations  6.3,  6.4,  6.6,  and  6.7, 

( 0-CHo ) 

II  3 

respectively 

and  predict  a and  x^3^  for  CH~OH  from  equations  6.5  and 

l|  ^ 

6.8.  The  results  of  this  fitting  along  with  the  experi- 
mental values  are  shown  in  table  9.  The  prediction  for 
a(CH3OH)  is  well  within  the  uncertainty  limits  expected 
for  these  data  (~  3-10$).  This  is  not  surprising,  since 
the  BAA  has  been  applied  successfully  to  a in  many 
cases. The  prediction  for  (CH^OH) , while  slightly 
outside  the  uncertainty  limits  of  the  experimental  value, 
it  is  still  not  unreasonable.  Thus,  it  appears  that 

bond  additivity  also  works  well  for  x^3^. 

II 
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The  situation  is  somewhat  more  complicated  for  the 
permanent  dipole  moment  (here  the  observable  p is  the 
magnitude  of  p,  the  dipole  moment  vector) . In  H^O  and 
(CH^)20,  symmetry  dictates  that  p must  lie  in  the  plane 
of  the  segments  and  along  the  bisector  of  the  angle  between 
them.  Thus  from  the  vectorial  addition  of  p,  we  find 


that 


p(O-H)  = 


h(h2o) 


p(o-ch3)  = 


h((ch3)2o) 


?COS  fro, 


(6.9) 


(6.10) 


where  0 is  the  angle ^ between  (0-H)  bonds  in  H20,  and  cp 
is  the  angle  between  (0-C)  bonds  in  (CH^^O.  These  bond 
moments  may  be  added  vectorially  at  the  proper  angle  to 
obtain  the  predicted  dipole  moment  for  CH30H.  The 
process,  angles  involved,  and  experimental  values  are 
shown  in  figure  7.  The  magnitude  of  the  predicted 
moment  is  within  10#  of  the  experimental  value,  although 
it  is  not  within  the  approximately  1#  experimental 
uncertainty.  In  addition,  the  orientation  of  the  molecular 
moment  is  within  10  degrees  of  the  experimentally 
determined  orientation.^ 

(2) 

In  applying  the  BAA  to  xj  for  these  molecules,  we 
are  hampered  by  the  lack  of  symmetry  about  the  plane 
bisecting  the  angle  between  the  two  segments  in  CH30H. 


□ 


Figure  7.  BAA  Applied  to  u for  H20,  (CH^O,  and  CH^OH. 
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!> 

From  equation  2.22,  we  see  th*fc--xj?i__is  the  z component 
U of  the  vector  xgT1T1»  where  the  z axis  is  along  the  dipole 

moment.  In  H^O  and  (CH-^O  this  vector,  which  we  shall 
call  Xy  , is  parallel  to  p.  by  symmetry,  so  for  these 
molecules  ^ 

u 


0 


(6.11) 


if 

a 


and  equations  analogous  to  6.9  and  6.10  hold: 


(h2o) 


2cos( 8/2) 


(6.12) 


x(2)(o-ch,) 
ii  3' 


X||2)((CH3)2°) 

2cos(cp/2) 


(6.13) 


ry  vector  addition,  we  can  determine  the  direction  and 
magnitude  of  x^.  '(CH^OH).  However,  the  lack  of  symmetry 
in  CH^OH  means  that  equation  6.11  does  not  hold. 

Instead,  we  must  use  the  more  general  relation 


(2)  - . y(2) 

' ' lei  K 


(6.14) 


where  the  direction  defined  by  h/|hI  is  that  determined 
from  experiment.  The  procedure  and  results  are  illustrated 
in  figure  8.  The  predicted  value  of  xjj2^  for  CH^OH  and 
the  experimental  value  agree  within  18$,  somewhat  outside 
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Figure  8.  BAA  Applied  to 
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the  experimental  uncertainty  of  about  6%. 

Thus,  we  have  seen  that  the  BAA  works  fairly  well 

in  correlating  electric  tensor  properties  of  this  simple 

set  of  molecules.  We  will  now  use  the  BAA  to  attempt  to 

correlate  the  properties  of  the  larger  set  of  halogenated 

methanes.  Those  molecules  of  this  set  which  possess 

permanent  dipole  moments  all  have  sufficient  symmetry 
(2) 

that  Xy  ' is  parallel  to  q,  providing  considerable 

simplification  in  the  calculation.  This  system  has 

three  types  of  bonds,  C-H,  C-Cl,  and  C-F,  and  thus  there 

are  three  bare  electric  tensors  of  each  order. 

Instead  of  fitting  to  molecular  observables  exactly  for 

several  molecules,  as  was  done  above,  we  have  chosen  bond 

2 

tensors  so  as  to  minimize  the  xv  (chi-squared  per  degree 
of  freedom)  defined  by 


(x(n)'_  x(n),2 
a 


(6.15) 


where  x^n^  are  the  experimentally  determined  observables 
i n 1 » 

values,  xv  ' the  fitted  values,  o are  the  experimental 

uncertainties,  v is  the  number  of  degrees  of  freedom,  and 

the  sum  is  taken  over  the  molecules  in  the  set.  This 

57 

standard  least-squares  procedure^  was  applied  to  a and 
X|j  ' of  the  halogenated  methanes,  and  the  results  are 
shown  in  figures  9 and  10,  respectively.  The  error  flags 
are  experimental  uncertainties.  As  can  be  seen,  the  fits 
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to  the  data  for  the  set  of  12  molecules  are  quite 

2 

good.  This  point  is  further  illustrated  by  the  Xv  of 
0.5  for  a and  17.4  for 

Application  of  the  BAA  to  p of  the  halogenated 
methanes  is  conceptually  no  different  from  the  previous 
case  and  yields  the  results  shown  in  figure  11  (this  set 
contains  only  9 molecules,  since  CF^,  CCl^,  and  CH^  have 
identically  zero  permanent  dipole  moments  by  symmetry). 

As  can  be  seen,  the  best  fit  to  the  data  is  poor.  The 
magnitudes  of  the  fitted  dipole  moments  differ  considerably 
from  the  experimental  values  (average  of  20  times 
experimental  uncertainty ), and  the  fit  also  fails  to 

reproduce  the  trends  of  p from  molecule  to  molecule. 

(2) 

XJj  ' can  be  fitted  in  exactly  the  same  manner  as  p,  and 
the  results  are  shown  in  figure  12.  Here  again,  the 
fitted  and  experimental  values  do  not  agree  well,  and  the 
molecule-to-molecule  trends  of  the  experimental  data  are 
not  reproduced  faithfully  by  the  fit. 

It  is  not  difficult  to  see  why  the  BAA  might  be 
expected  to  fail  in  a molecular  series.  Consider  two 
polarizable  bodies  A and  B separated  by  a distance  R as 
in  figure  13,  having  permanent  dipole  moments  pA  and  pfi 
of  about  the  same  magnitude,  and  isotropic  linear 
polarizabilities  aA  and  aB,  also  comparable  in  size.  The 
electric  field  at  A due  to  the  dipole  moment  of  B will  be 
of  the  order  p^R-^,  and  therefore  the  dipole  moment  ApA 
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Figure  13.  Interaction  Between  Neighboring 
Polarizable  Dipolar  Bodies. 
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_ O 

induced  in  A due  to  will  be  approximately  p^o^R  . 

-24  8 

Molecular  polarizabilities  are  of  the  order  10  cm  , 

-8 

and  typical  inter-segment  distances  in  molecules  are  10  cm, 

_3 

so  the  factor  aAR  J is  of  order  1,  and  AnA  is  of  the  order 

of  a segment  dipole  moment.  Clearly  the  same  sort  of 

argument  applies  to  other  orders  of  segment  tensors.  In 

the  BAA,  the  anisotropy  of  these  intersegment  interactions, 

and  their  variation  from  molecule  to  molecule,  are 

ignored.  Indeed,  it  is  more  surprising  that  BAA  should 

work  so  well  for  a and  than  that  it  should  fail  for 

(2)  11 

u.  and  xv  . In  addition,  BAA  has  been  quite  successful 
11 

in  correlating  hyperpolarizabilities  in  crystals-*  and 
■30 

liquids.  It  has  been  argued  that  cancellation  among 

inductive  terms,  due  to  statistical  isotropy  in  liquids 

or  symmetrical  arrangement  of  neighboring  segments  in 

crystals,  might  be  responsible  for  the  success  of  BAA  in 

these  media.  Such  isotropy  and/or  symmetry  does  not, 

however,  apply  to  a segment  in  an  isolated  molecule. 

84 

Attempts  by  some  early  authors-^  to  account 
approximately  for  inductive  effects  among  the  C-H,  C-Cl, 
and  C-F  bond  dipole  moments  of  the  halogenated  methanes 
met  with  limited  success  in  predicting  the  molecular 
dipole  moments.  However,  iterative  methods  were  used 
in  that  work.  We  will  consider  a self-consistent  approach 
to  calculating  inter-segment  interactions  in  molecules. 
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6.2  The  Interacting  Segment  Model 


We  wish  to  develop  an  Interacting  Segment  Model  (ISM) 

which  considers  bond  interaction  contributions  to  the 

electric  tensor  properties  of  structurally  similar 

molecules,  without  completely  discarding  the  simplicity 

of  the  BAA.  A molecular  group-dipole  interaction  model 

bearing  formal  similarity  to  our  ISM  has  recently  been 

■zq 

developed  by  Sundberg.  ^ 

All  segment  electric  tensor  properties  in  the  BAA 
have  two  characteristics  in  common.  These  are  additivity, 
as  expressed  by  equation  6.1,  and  transferability,  meaning 
that  bond  properties  are  independent  of  molecular 
environment.  In  the  ISM,  we  assign  these  characteristics 
to  two  different  sets  of  tensors.  We  suppose  that  in 
the  absence  of  intersegment  interactions,  a segment  i 
can  be  characterized  by  a set  of  "bare"  segment  parameters 
J±j_,  these  parameters  are  to  be  thought 

of  as  the  electric  tensor  properties  of  an  isolated 
segment  so  that,  for  example,  a C-X  bond,  X = halogen,  has 
rotational  symmetry  about  the  bond  axis.  These  bare 
parameters  are  defined  to  be  transferable  from  molecule 
to  molecule.  When  intersegment  interactions  are  included, 
the  segment  i is  characterized  by  a set  of  "dressed" 
parameters  a^,  these  tensor  properties 

have  been  modified  by  interactions  with  neighboring 
segments,  and  will  possess  a symmetry  dictated  by  the 
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arrangement  of  surrounding  segments,  which  will  in  general 
be  lower  than  that  of  the  bare  parameters. 

These  tensors  are  assigned  the  characteristic  of  additivity 
and  so  conform  to  equation  6.1,  which  now  takes  the  form 


where  the  sum  is  over  all  segments  in  the  molecule. 

It  is  appropriate  at  this  point  to  adopt  the 

dispersionless  limit,  that  is,  to  asstime  that  all 

are  equal  to  their  dc  values.  As  we  have  seen,  the  fact 
(2)  (q ) 

that  Rv  ' and  RVJy  are  very  nearly  3.0  implies  that 
dispersion  is  not  important  in  these  molecules,  and 
by  neglecting  its  effect,  we  will  realize  a significant 
simplification  in  the  mathematics.  In  this  limit,  we 
can  write  the  dipole  moment  of  segment  i in  the  presence 
of  an  external  electric  field  E as 


Ui(E)  = Hi  + “iH  + ii2^2  + £j3^E3+  ... 

= S Ki11^  (6.17) 

where  En  is  a tensor  of  rank  n,  the  direct  product  of  E 
with  itself  n times,  and  is  the  inner  product 

of  the  (n+l)th  rank  tensor  with  En.  Thus,  equation 

6.17  is  equivalent  to  equation  2.1  for  dc  field,  except 
that  it  is  applied  here  to  a segment,  rather  than  to  a 
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complete  molecule.  One  can  also  express  ^(E)  in  terms 
of  the  total  electric  field  at  segment  i,  consisting 
of  the  external  field  E and  the  fields  at  segment  i due 
to  permanent  and  induced  dipoles  at  neighboring  segments. 
Since  the  interactions  are  included  in  the  electric  field, 
this  expression  involves  bare  segment  tensors  and  is 
written  as 

4^  (E)  = y.i  + a1Ei  + x{2)E^  + Xp^E^  + ... 

(6.18) 

In  equations  6.17  and  6.l8,  we  have  assumed  that  all 
intersegment  interactions  are  electrostatic  in  origin, 
and  we  have  ignored  all  non-classical  short-range  inter- 
actions. To  obtain  simple  expressions  for  E^,  we  further 
assume  that  all  tensors  for  a given  segment  are  localized 
within  the  segment  at  a single  point,  called  the  "active 
center",  and  that  contributions  from  segment  electric 
multipoles  of  order  higher  than  dipole  can  be  ignored. 

This  point  localization,  which  was  not  required  by  the 
BAA,  finds  some  justification  in  localized  molecular 
orbital  calculations,  but  it  is  still  a major  simplifi- 
cation. Such  drastic  approximations  can  only  be 
justified  by  the  ability  of  the  model  to  reproduce 
experimentally  observed  molecular  tensor  properties,  given 


= S x[n)E? 
n=0 
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a single  set  of  bare  parameters  which  are  transferable 
from  molecule  to  molecule.  With  these  approximations, 
the  electric  field  AE^1^  at  segment  i due  to  the  dipole 
Hi  (ill)  at  segment  j is  T^jijCE)  or,  more  explicitly. 

Hereinafter,  Latin  indices  refer  to  segment  numbers  and 
Greek  indices  stand  for  spatial  coordinates,  and  summation 
is  implied  over  repeated  Greek  indices  only.  The  tensor 
is  defined  by 

Tijtr,  * - fienRtj2]  Rij5’  3/1 

= 0 , j=i  (6.20) 

where  R,  . is  the  distance  vector  between  active  centers 
i and  j,  and  6.  is  the  Kronecker  delta. 

The  total  electric  field  at  segment  i can  then  be 
written 


Ei  - E + E T1jEj(E) 


(6.21) 


Substituting  equation  6.17  into  equation  6.21  yields  an 
expression  for  in  terms  of  the  applied  field  E,  the 
dressed  tensors  and  the  geometrical  factors  T^. 

Further  substitution  of  equation  6.21  into  equation  6.18 
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results  in  a polynomial  in  E which  can  be  compared  term- 
by-term  with  equation  6.17  to  yield  expressions  for 
in  terms  of  the  T^,  of  a11  orders  p,  and  £(in)  for 

n < m.  Before  writing  down  these  results,  we  introduce 
further  notation  conventions: 

1.  The  expression  <T*^n)\  stands  for  the  sum 
over  j of  the  inner  product  of  and 

2.  All  products  of  such  terms  (including  the  m-fold 
product  ( T*^n^ >^)  are  direct  products; 

3.  The  product  of  a bare  parameter  with  a term 
such  as  one  of  those  in  1 and  2 above  is  an 
inner  product. 

Using  these  conventions,  results  for  the  dressed  parameters 
are  written  as 

i±i=  Jii+  x[2^<T*m.>2  + + ...,  (6.22) 

«I=  (0^+  2x{2)<T-ji>i+  3x(3)<T4i>2+  . . . ) (l+< T*&1  i) , 

(6.23) 

ii2)=  “i<ri(2)>i+  x[2)C2<T-M;1><T-^2^i+(l+<T.a>i)2] 

+ 3^[3)<T4>1C<T-H>i<T*^2^>i+(1+<T4>i)2;i+  • ••*  I] 

(6.24)  | 


• . . . 

. ,•  v • 

— ...  . — ■ — —■ 


i 


ii3)=  £i<Ti(3)>i+  2x[?)[<T-4>i<T-^3)>  + (l+<T*a>i)<T'i(2)>i] 
+ 43)C3<T4>i<i,^3H-t6<T4>i(l+<T*a>i)<T'i^Hi 


+ ( l+< T* i ) 3]  + ...  . 


(6.25) 


Terms  not  explicitly  included  in  equations  6.22-6.25  are 
those  containing  with  n > 3.  The  dimensional 

arguments  in  the  previous  section  and  in  section  2.1 
would  seem  to  suggest  that  for  all  orders  n there  will  be 
terms  in  which  will  be  comparable  in  magnitude  with 

the  leading  terms  in  these  equations.  However,  we  adopt 
the  hypothesis  that  higher-order  are  overestimated  by 

the  dimensional  arguments.  Also,  cancellation  in  the 
tensor  sum  over  bonds  <T*  where  is  an  even- 

order  electric-tensor  which  vanishes  in  a centro- symmetric 
system,  will  reinforce  this  tendency  (contributions  to 

which  contain  x/m^»  m > n,  must  have  at  least  a 
factor  (T-J.)^'11). 

The  task  that  remains  is  to  determine  sets  of  dressed 

and  undressed  molecular  tensors  from  molecular  observables, 

using  the  expressions  developed  in  this  section.  We  now 

consider  this  problem  for  the  permanent  dipole  moments  ^j., 

linear  polarizabilities  a,  and  second-order  polarizabili- 
(2) 

ties  ’ of  the  halogenated  methanes. 


6.3  Application  of  the  ISM  to  Electric  Tensor  Properties 


of  the  Halogenated  Methanes 


Equations  6.22-6.25  form  a coupled  non-linear  system. 
Our  goal  is  to  find  a set  of  bare  bond  tensors  which,  when 
substituted  into  these  equations,  reproduce  the  molecular 
observables.  We  will,  however,  make  some  simplifying 
assumptions  before  doing  so. 

Our  first  step  is  to  exclude  x^)  from  the  fitting 
procedure,  eliminating  equation  6.25  from  the  set.  Since 
only  the  isotropic  parts  of  the  molecular  x^^s  are  known, 
it  is  felt  that  anisotropies  of  the  bare  parameters  x|^ 
would  be  ill -determined  by  fitting.  Thus  the  bare  x[^ 
will  be  considered  to  be  isotropic,  with 


4zl2<c-x>  - • x - f’ci’h 


(6.26) 


Secondly,  we  invoke  our  hypothesis  that  terms  con- 
taining tensor  sums  <T*^^n)>  will  be  reduced  in  magnitude 
owing  to  cancellation  in  the  sums.  We  will  ignore  any 
terms  containing  even-order  polarizabilities  in  expressions 
for  with  odd  n,  and  we  will  discard  terms  containing 

more  than  one  even-order  factor  in  expressions  for 
with  even  m.  Thus,  the  equations  for  will  contain  no 
even-order  factors,  and  terms  in  the  equations  for  ^ and 
£(2)  will  have  no  more  than  one  even-order  factor.  With 
this  restriction,  we  rewrite  equations  6.22-6.24: 
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£i  = Hi  + SLi<T4>i  i (6.27) 

-i  = -i  + > (6.28) 

x[2)  = X^2)  + £i<T*i(2)>i+  x[2)<T-a>1(2+<T4>1) 

+ 3x[3^<T4>i(l+<T-ft>1)2  . (6.29) 

In  employing  this  approximation,  we  have  gained  a 

significant  advantage.  While  equations  6.22-6.24  were 

coupled,  6.27-6.29  are  independent  in  the  following  sense: 

can  be  selected  to  fit  molecular  a values  using 

equation  6.28,  without  reference  to  ji  or  ^ can  then 

be  found  so  as  to  fit  molecular  values  by  means  of  equation 

6.27,  using  the  calculated  in  the  previous  step}  and 
(o)  ( 2 ) 

finally,  ^ ' can  be  selected  to  fit  molecular  xjj  values, 
making  use  c"  equation  6.29  and  the  jl^,  a^,  and  a.  already 
determined.  The  steps  in  this  procedure  are  elaborated 
upon  below. 

Before  equations  6. 27-6.29  can  be  used  for  calcula- 
tions, the  T^j  must  be  evaluated,  i.e.,  active  center 
locations  must  be  selected.  We  have  arbitrarily  chosen 
as  the  distance  from  the  carbon  nucleus  to  the  C-X  bond 

-3K 

active  center  the  average  C-X  bond  length,  ^ and  for 
angles  we  use  experimentally  determined  bond  angles. 

These  values  are  shown  in  table  10.  We  have  explored 


Table  10. 


Molecule 

Bond 

Angles  a 

CCl^ 

1C1-C-C1 

= 

109°28' 

b 

CHC13 

Z.H-C-C1 

= 

1080 

CH2C12 

1H-C-H 

= 

112° 

Z.C1-C-C1 

= 

111°48' 

CH3C1 

tCl-C-H 

= 

108° 

ch4 

Z.H-C-H 

= 

109°28' 

b 

CH3F 

1 

0 

1 

= 

108o49’ 

CH2F2 

Z.H-C-H 

= 

111°54’ 

1F-C-F 

= 

108°30’ 

CHF3 

Z.H-C-F 

=r 

110°8' 

CFi/ 

Is- 

1 

0 

1 

= 

109°28' 

b 

CC1F3 

Z.C1-C-F 

= 

110°20' 

cci2f2 

^Cl-C-Cl 

= 

108°30' 

Z.F-C-F 

= 

109°30* 

CC13F 

LF -C-Cl 

= 

107°12 ' 

Bond 

C-F 
C-Cl 
C-H 

aFrom  reference  35. 
^Tetrahedral  symmetry. 
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other  reasonable  choices  for  bond  radii,  with  results 
qualitatively  similar  to  those  shown  here. 

It  is  desirable  to  rewrite  equation  6.28  so  that 
is  not  functionally  dependent  upon  a^.  We  write 

“ilijlij  - Si  (6.30) 

where  is  an  inner  product  and  the  are  tensors 

such  that 


= 6446- 


(6.31) 


We  then  take  the  inverse  M of  the  matrix  quantity  {...} 
in  equation  6. 30  to  obtain  the  expression 


!i  - a^jSj  (6.32) 

This  is  identical  to  the  result  of  Applequist  _et  al. 

obtained  in  investigations  of  the  polarizability  of 

methane.  Their  B^.  correspond  to  the  negative  of  the 

inner  product  M^a^.  Equation  6.16  is  used  to  obtain 

tensors  a' , and  standard  non-linear  least-squares  fitting 
37 

techniques-^  are  used  to  adjust  to  get  the  best  fit 
of  the  a*  to  molecular  observables.  In  addition  to  the 
isotropic  polarizabilities  a defined  by  equation  6.2, 
there  are  two  other  observables  derivable  from  a.  These 


|p> 
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p 

are  h , defined  by 


^ = I^t(a«-0;yy)2+(ayy'a^)2+(a^'c‘xx)2:l'  <6-33> 


obtained  from  measurements  of  the  depolarization  of 
Rayleigh  scattering,  and  a quantity  which  we  shall  call 
AaK,  defined  by 


4aK 


(6.34) 


which  is  extracted  from  measurements  of  the  temperature 
dependence  of  the  dc  Kerr  effect.  For  this  last  observ- 
able, the  z axis  is  along  the  dipole  moment,  so  that  in 
order  to  obtain  corresponding  Aa^  values  from  the  a',  it 
is  necessary  to  relate  this  axis  to  molecular  geometry. 
Those  halogenated  methanes  of  this  set  which  possess 
dipole  moments  have  either  C2v  or  C^v  symmetry,  and  the 
dipole  moment  must  lie  along  the  major  symmetry  axis. 

Once  the  are  determined,  the  are  fixed  and 
equation  6.27  may  be  manipulated  to  yield 

• (6.35) 

<j 

A simple  linear  least-squares  procedure  is  then  used  to 
fit  the  resulting  p.'  to  the  molecular  p.  values.  In  order 
to  fit  p.*  to  p.  unambiguously,  it  is  necessary  to  know  the 


9b 


sense  of  the  dipole  moment  in  relation  to  molecular 

geometry.  It  is  taken  here  to  be  in  the  sense 

~FnCCl^_n,  +HnCCl^_n,  or  +HnCF4_n*  which  is  consistent 

4l  + 

with  the  experimentally  determined  polarities  HCF^ 
and  +H3CC1“. 

The  thus  determined,  along  with  the  tensors 
and  the  previously  determined  a^,  appear  in  the  expression 
obtained  by  solving  equation  6.29  for 

x[2)  = E Mij(*j2)+  ^j3)<l4>j)(1+<T4>j)2-  (6.36) 

Again,  all  that  is  required  is  a linear  least-squares  fit 

(2)  (2) 

of  the  resulting  x ^ bo  the  molecular  xjj  ' values. 

The  results  of  fitting  these  functions  to  experimental 

data  are  shown  in  figures  14,  15,  and  16,  which  may  be 

compared  with  the  BAA  fits  in  figures  11,  9>  and  12, 

respectively.  The  qualitative  nature  of  the  fit  to  the 

dipole  moments,  shown  in  figure  14,  is  not  dramatically 

better  than  that  provided  by  the  BAA.  There  is,  however, 

some  improvement,  particularly  for  the  small  dipole 

moments  of  the  CC1  F>.  „ molecules.  The  fit  to  the 

isotropic  part  of  a,  shown  in  figure  15,  remains  excellent. 

In  addition,  there  is  a qualitatively  good  fit  to  most  of 

the  A a and  h2  data.  The  x^  fit  shown  in.  figure  16  is, 

R II 

like  the  4 fit,  not  strikingly  improved  over  bond 
additivity.  Again,  however,  there  is  some  degree  of 
qualitative  improvement,  this  time  in  the  CHnCl^_n  and 


Figure  16.  ISM  Applied  to  x ' for  the  Halogenated  Methanes 


Terms  contributing  to  fit 


Values  of  Aa^.  from  the  various  terms  are  shown  in 

table  13.  The  and  ^(T-a^  terms  both  make  large 

contributions  to  this  observable.  The  contributions 
(2) 

from  the  ^ ' term,  which  was  not  included  in  the  fit, 
are  large,  particularly  for  the  CHnF^_n  molecules.  The 
contributions  from  the  term  involving  are  much 

smaller. 

(2\ 

Contributions  to  the  xv  fit  from  terms  in  the 

11 

expression  for  are  shown  in  Table  14.  Contributions 

are  generally  huge  compared  to  the  X^'  values.  In 

addition,  the  neglected  terms  involving  x[2^  are  large. 

The  terms  involving  x^)  which  were  not  included  in  the 

fit  are  small  except  in  the  case  of  CH2C12.  All  of  the 

contributions  are  large  for  this  molecule.  This  situation 

is  not  surprising  when  we  look  at  the  x^2^  for  the 

individual  bonds.  The  dressed  bond  tensors  are  large,  but 

they  cancel  to  less  than  a part  in  200  when  combined  to 

(2) 

form  a molecular  xjj  . However,  the  individual  terms  do 
not  each  cancel,  and  so  appear  as  large,  conflicting 


n 

0 


molecular  contributions.  In  the  case  of  CH2C12,  for 

(2) 

example,  the  neglected  contribution  involving  ^ ' is 

less  than  10$  of  X.[2^  in  each  of  the  two  bond  types  in 

the  molecule,  yet  when  summed  over  the  molecule,  it  is 

over  6 times  the  total  x^2^'.  This  situation  occurs,  to 

II 

some  degree,  in  all  of  the  molecules  in  this  set,  not  only 
for  x^2^  but  for  4 and  AaR  as  well.  In  view  of  the 


Table  13. 
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Definition  of  xyc ' Terms  for  Table  14 


b)  x[2)<t4>1 

c)  a^T'x/2^ 

d)  3x[3)<T-4>i(l+<T-a>i)2 

e)  2xi[2)<T*U>i<l-X(2)>i 

f)  3x[3)<t4>2<t-x(2)>1 
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degree  of  cancellation,  it  is  somewhat  surprising  that 
the  fits  are  as  good  as  they  are. 

In  summary,  we  have  constructed  a model  of  electric 
tensor  properties  in  which  inter-bond  interactions  are 
taken  into  account.  We  include  only  electrostatic  dipolar 
interactions,  ignoring  short-range  non-classical  forces 
and  all  higher  order  multipole  moments.  Further,  we 
assume  that  all  tensor  properties  of  a bond  are  located 
at  a single  point,  the  position  of  which  is  arbitrary, 
although  consistent  with  physical  intuition.  The  model 
meets  with  several  difficulties  in  implementation.  Some 
of  the  observables  are  small  differences  resulting 
from  incomplete  cancellation  of  large  bond  tensors.  This 

is  the  case  for  4 and  and,  to  some  degree,  for  AaK 

0 " / p ) 

and  h . In  addition,  the  x ^ and  higher  order  quantities 
are  increasingly  complicated  and  difficult  to  calculate. 

To  eliminate  coupling  between  the  coefficients,  it  was 
necessary  to  truncate  the  expressions  for  the  dressed 
tensors  in  a consistent  manner,  leaving  what  is  in  some 
sense  a "first  level"  of  interactions. 

Despite  the  foregoing  difficulties,  our  model  still 
provides  a fit  to  the  experimental  data  which  is  somewhat 
better  than  that  of  the  BAA.  In  particular  a,  which  BAA 
fits  very  well,  is  fit  equally  well  by  the  ISM.  Thus 
we  have  established  that  the  fact  that  isotropic 
properties  are  well  fit  by  the  BAA  is  not  sufficient  to 
imply  that  bond-bond  interactions  have  a negligible 


effect  on  molecular  electric  tensor  properties.  A 
posteriori  evaluation  of  neglected  terms  indicates 
some  trend  toward  decrease  in  the  "odd-parity"  contribu- 
tions involving  and  such  terms  involving  are 

shown  to  be  much  less  important. 

In  view  of  these  encouraging  signs,  we  plan  to 
extend  our  calculations  to  include  the  terms  neglected 
here,  including  equation  6.25.  We  feel  that  there  is 
reason  to  hope  that  such  a fitting  procedure,  although 
formidable  to  implement,  would  arrive  at  a self-consistent 
set  of  bare  bond  parameters  yielding  a good  fit  to  the 
experimental  data  and  giving  a satisfactory  picture  of 


inter-bond  interactions. 


CHAPTER  7 


u 


m 

a 
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SUMMARY  AND  CONCLUSIONS 

The  objectives  of  the  research  presented  here  have 

been:  measurement  of  the  third  order  polarizability  x^3^ 

and,  where  appropriate,  the  second  order  polarizability 
(2) 

Xv  ' of  a number  of  simple  molecules  of  potential  or 
demonstrated  theoretical  and  experimental  interest,  using 
the  technique  of  dc-electric-f ield-induced  optical 

(2) 

second-harmonic  generation  (dcSHG);  measurements  of  x ' 
and  X(3)  for  those  molecules  of  the  set  CX^Y^  n*  n = 1-4, 
X,Y  = F,C1,H,  for  which  measurements  did  not  already 
exist;  and  development  of  a model  to  describe  trends  in 
the  dipole  moment,  linear  polarizability,  and  hyper- 
polarizabilities for  series  of  structurally  related 
molecules,  including  the  effects  of  inter-bond  inter- 
actions, with  application  to  the  halogenated  methanes. 

A review  of  second-  and  third-order  electric  tensor 
properties  was  presented.  Also  reviewed  was  the  dcSHG 
process  in  gases  for  a particular  electric  field 
configuration,  with  the  treatment  extended  to  include 
anomalously  dispersive  gases. 

The  measurements  were  performed  using  an  apparatus 
consisting  of  a pulsed,  Q-switched  ruby  laser  and 
associated  optics  and  electronics,  along  with  an 
experimental  cell  and  gas  handling  system  designed  to 
permit  temperature-dependent  measurements,  in  the  gas  phase 
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on  substances  which  are  liquids  at  room  temperature. 

The  hyperpolarizabilities  measured  for  the  assorted 
small  molecules  were  compared  with  other  gas-phase 
experimental  data  and  with  theoretical  results,  wherever 
such  data  were  available.  Our  values  agreed  acceptably 
with  most  other  experimental  data.  In  cases  where  large 
discrepancies  occurred,  it  was  felt  that,  because  of  the 
nature  of  the  experiments,  our  data  were  the  more 

(2) 

reliable.  Agreement  with  theoretical  estimates  for  xv  ’ 
was  poor,  but  calculations  by  various  authors  also  dis- 
agreed strongly  with  each  other.  A quantum-mechanical 
approximation,  formulated  by  Dawes  to  predict  ' for 
the  inert  gases,  was  applied  to  N^,  and  O^,  but 

(2) 

agreement  with  experiment  was  poor.  Measurements  of  xv 

and  X(3)  were  also  made  for  several  halogenated  methanes, 

and  these  results,  too,  were  compared  with  other 

experimental  data.  Agreement  with  other  x^3^  data  is 

(2  ^ 

excellent.  The  xv  from  other  gas-phase  measurements 
tend  to  be  larger  than  our  values,  but  for  the  most 

part  there  is  no  serious  inconsistency.  Agreement  with 

( ?\ 

\x  > values  extracted  from  liquid  phase  measurements  is 
poor,  but  this  is  not  unexpected,  since  the  derivation 
of  values  from  liquid  measurements  is  indirect.  There 
were  no  theoretical  x^2^  or  x^3^  values  available  for 
these  molecules.  However,  it  is  anticipated  that  new 
theoretical  techniques  will  lead  to  accurate  calculations 


■ 
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for  many  of  the  molecules  for  which  we  have  made 
measurements.  Our  data  will  provide  a stringent  test  of 
the  validity  of  such  calculations. 

Correlation  of  molecular  electric  tensor  properties 

among  molecules  of  similar  structure  was  discussed.  The 

Bond  Additivity  Approximation,  or  BAA,  was  applied  to 

4,  a,  and  of  the  set  HpO  , (CH^)pO,  and 

CH^OH,  and  was  found  to  agree  well  with  experimental  data. 

This  was  also  the  case  when  the  BAA  was  applied  to  a and 

X(3)  for  the  halogenated  methanes.  However,  the  BAA  did 

not  reproduce  well  the  variation  from  molecule  to  molecule 

(2) 

of  the  experimental  p.  or  x v values.  A model  of 
molecular  electric  tensor  properties  including  the  effects 
of  dipolar  electrostatic  interaction  among  bonds  was 
introduced.  A simplified  version  of  the  model  was 
applied  to  a,  p.,  and  xjj^  of  the  halogenated  methanes, 
with  a modest  improvement  over  BAA,  including  the  case  of 
a,  in  which  BAA  worked  well.  We  conclude  that  the  success 
of  BAA  in  dealing  with  isotropic  properties  does  not 
imply  that  interactions  among  bonds  can  be  neglected.  It 
was  found  that  higher  order  terms  in  our  model  tended 
to  make  smaller  contributions  to  the  observables.  It  is 
felt  that  the  results  warrant  the  effort  to  implement  the 
full  bond-interaction  model,  including  all  relevant 


) 


I 


contributions. 
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APPENDIX 

DCSHG  IN  ANOMALOUSLY  DISPERSIVE  GASES 
The  development  in  chapter  two  of  an  expression  for 
second  harmonic  generation  as  a function  of  experimental 
geometry  required  a choice  of  sign  to  be  made  for  the 
wave-vector  mismatch  Ak.  The  choice  leading  to  equation 
2.28  is  Ak  < 0,  which  is  appropriate  for  gases  of  normal 
dispersion.  We  demonstrate  here  that  the  occurrance  of  a 
positive  Ak  would  be  readily  detectable  in  a binary 
mixture  experiment  with  a normally  dispersive  gas. 

The  wave-vector  mismatch  Ak  enters  into  the  derivation 
of  equation  2.28  in  the  integral 


I 


iAkz' 

e 


1 


o'\2 


dz ' 


(Al) 


where  the  limits  of  integration  are  appropriate  to  harmonic 
generation  in  an  infinite  medium,  which  is  a good  approxi- 
mation here.  For  Ak  < 0,  the  integration  is  performed  by 
the  method  of  residues,  considering  as  the  contour  the 
real  axis  and  a semicircle  at  infinity  in  the  lower  half- 
plane. The  only  contribution  to  this  integral  is  the 
single  pole  in  the  lower  half-plane  due  to  the  dc  field. 

For  Ak  > 0,  the  contour  of  the  integral  must  be  completed 
in  the  upper  half-plane,  in  which  both  the  dc  field  and 
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the  optical  field  have  poles.  This  integration  can  be 
done,  yielding  an  expression  for  S>,  the  value  of 
E2aV(F^)2  for  Ak  > 0,  which  is 


= S 


* = < (Zo  + (V)y)  C 


1 ;d2e-(b-d)|Ak| 


2de"(b"d) lAkl/2[z  sin(|k|z  ) 


- (b^)cos(|Ak|zo)]  + z2  + (^)2} 


(A2) 


where  S<  is  the  expression  for  E2cD/(lP)2  from  equation 
2.28.  It  should  be  noted  that  only  the  magnitude  of  Ak 
appears  in  equation  2.28.  The  expression  for  S>  has 
independent  maxima  at  zq  = 0 and  b = d,  as  does  S<,  so 
again  it  is  appropriate  to  work  at  these  values  of  b and 
zq.  It  would  appear  at  first  glance  that  equation  A2  has 
a singularity  at  zq  = 0 and  b = d.  However,  with 
zq  = 0 approaches  a finite  limit  as  b approaches  d,  and 
this  is  the  same  limit  approached  by  with  b = d as 
z^  approaches  zero.  The  limit  is 

S>  = (1  + d | Ak | )2S<  . (A3) 

{ 


With  this  expression,  we  are  now  in  a position  to  predict 
the  signal  for  binary  mixtures  of  gases. 
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Suppose  that  we  have  two  gases,  labeled  1 and  2, 
such  that  (AkQ)2  < 0,  and  the  sign  of  x^2  (which  we  shall 
call  simply  x2)  is  known.  We  will  refer  to  the  |x^| 
measured  for  gas  1 by  the  method  described  in  section  3.2a 


as  1x^1  . If  (Ak^  is  negative 
Ix-J  = Ix-l  I , (AkQ)i  < 0 


However,  if  (AkQ)1  is  positive,  then 


Ix-J  = I Xx  | /3  , (Ako)1  > 0 


The  factor  of  three  arises  from  the  fact  that  the  measure- 
ment is  done  with  p = pQ.  The  optimum  density  oq  is 
defined  by 


d|AkQ| 


where 


|Akol  = | Ak | /p  . 


From  equation  A3  we  see  that  for  positive  AkQ,  measurement 
at  pQ  leads  to  a factor  of  9 greater  signal  for  a given 
|X®I  than  in  the  case  of  negative  AkQ.  This  factor  is 
not  accounted  for  in  the  procedure  of  section  3.2a. 


> • K~r  , -v  • , • 
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For  a mixture  of  gases,  the  effective  third-order 
polarizability  is  given  by 


XM  “ 


X1P1  + X2P2 
pl  + P2 


and  the  wave-vector  mismatch  AkM  by 


' pl<iko>l  + P2(lko>2 


For  AkM  < 0,  equation  2.28  can  be  used  to  find  the  form  of 
the  signal  S,_,  generated  in  the  mixture: 

M< 


SM<  = CIX1P1  + X2P2  * 2 exP(*dlAkM^  ‘ 


where  C is  a constant  independent  of  the  gases. 


(A10) 


In  practice,  a fixed  amount  of  gas  2 is  placed 
in  the  cell,  and  the  variation  of  the  signal  SM  is 
measured  as  the  density  p.^  of  gas  1 is  increased. 
Combining  equations  A3,  A4,  A5 , and  A10  we  get 


SM  = So  ,X1P1+  X2P2|2exp(-d|AkM|)  ' (Akoh<0 
SM  = So^xlpl+  X2P2|2exp(-d|AkM|)  , (Ak^^O. 


^ (Ako>2 
Pi  — — P2 

1 (AkoU  2 
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SM  = Sol^lpl+  X2p2|2exp(-dkkMl  Hi+dlAknl  ),  (Ako)1>0, 

( Ako ) 2 


pl> 


(Ako\ 


(All) 

where  SQ  is  a factor  normalizing  the  computed  SM  values 
to  the  measured  signal  value  at  p^=  0.  SM  can  be  plotted 
for  four  different  cases,  as  follows: 

1.  (AkQ)1<0,  and  x2  of  opposite  sign.  As  p-j^ 

increases,  the  signal  will  decrease  until 
cancellation  occurs  at 


pl  = 


(A12) 


The  signal  will  increase  again,  and  then 
decay  exponentially. 

2.  ( AkQ ) -|< 0 , x^  and  XP  same  sign.  As  p-j^  increases, 
signal  increases  to  a maximum  and  then  falls 
off  exponentially.  We  choose  p^  to  maximize 
the  difference  in  signal  between  case  1 and 
case  2 at  the  expected  cancellation  density. 

As  Finn^  has  shown,  this  density  is 


|fe|  |(lk0)1|  + |(4k0)2 


(A13) 
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3.  (Ako).j>0,  x1  and  x2  of  opposite  sign.  As  p1 

increases,  the  signal  may  increase  at  first, 
but  will  then  decrease  to  zero  at 


(A14) 


I 

E 

0 


It  then  increases  again,  and  finally  falls  off 
exponentially.  Note  that  the  expected  cancellation 
density  is  three  times  that  of  case  1. 

4.  (Ako)1>0,  x1  and  x2  of  the  same  sign.  As  p1 

increases,  the  signal  increases  rapidly,  and 
reaches  a maximum  much  larger  than  that  of  case  2. 
Beyond  the  maximum,  the  signal  decreases 
exponentially. 

Plots  of  the  expected  signal  for  these  four  cases, 
along  with  experimental  data,  are  shown  in  figure  17.  Gas 
2 in  both  plots  is  methane.  In  the  upper  plot,  gas  1 is  water 
vapor,  and  in  the  lower  plot,  gas  1 is  methylene  chloride 
vapor.  The  data  show  unambiguously  that  x^(H20)  is  of 
opposite  sign  to  x^(CH^),  x^(CH2C12)  is  of  the  same  sign 
as  x^(CH^),  and  AkQ  for  both  gases  is  negative. 

We  have  shown  that  a binary  mixture  experiment  can  be 
used  to  determine  the  sign  of  A kQ  for  a gas  simultaneously 
with  the  determination  of  the  sign  of  x®  • Such 
experiments  have  shown  that  AkQ<0  for  all  of  the  gases 
that  we  have  investigated. 
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